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ELECTRON  LEAKAGE  CAUSED  BY  STOCHASTIC  ELECTRON  MOTION 
IN  MAGNETICALLY  INSULATED  DIODES 


M.  P.  Desjarlau  and  R.  N.  Sudan 


Laboratory  of  Plasma  Studies 
Cornell  University 
Ithaca,  New  York  14853 


Abstract 

The  diffusion  of  electrons  and  the  associated  leakage 
current  is  investigated  for  magnetically  insulated  sys¬ 
tems.  The  diffusion  results  from  stochastic  motion  of 
the  electrons  in  a  weakly  nonuniforin  electric  field.  The 
model  system  is  planar  and  periodic  in  the  electron  drift 
direction.  Leakage  currents,  normalized  to  the  Child- 
Langtnuir  current,  are  obtained  as  a  function  of  the  in¬ 
sulation  parameter  B/Br.  Very  good  agreement  with 
experimental  results  is  obtained  for  electric  field  pertur¬ 
bations  on  the  order  of  5  to  7%  of  the  uniform  equilib¬ 
rium  field. 


Introduction 

Magnetic  insulation  of  electrons  in  crossed-ficld  devices 
plays  a  crucial  role  in  the  generation  of  high-power  ion 
beams  and  the  operation  of  magnetically  insulated  trans¬ 
mission  lines.  Successful  application  of  these  devices  re¬ 
quires  that  electron  losses  be  kept  to  tolerable  levels. 
Simple  models  for  the  electron  sheath  that  assume  sym¬ 
metry  in  the  drift  direction  predict  complete  insulation 
for  sufficiently  large  magnetic  fields.  Breaking  of  this 
symmetry,  whether  by  imperfections  of  construction  or 
the  spontaneous  growth  of  fluctuations,  may  result  in 
stochastic  wandering  of  the  sheath  electrons  and  the  gen¬ 
eration  of  an  associated  leakage  current. 

In  an  earlier  work  |lj,  we  investigated  the  transition 
to  stochasticity  for  electron  motion  in  the  field  configura¬ 
tion  of  a  magnetically  insulated  diode;  action-angle  vari¬ 
ables  were  introduced  to  facilitate  the  nonlinear  analysis 
and  the  criterion  for  stochasticity  was  established.  In 
this  paper,  we  calculate  the  resulting  diffusion  and  leak¬ 
age  currents.  Using  quasilinear  theory,  the  diffusion  co¬ 


efficient  is  calculated  for  a  diffusion  equation  relevant  to 
electrons  in  the  globally  stochastic  regions.  The  steady- 
state  diffusion  equation  is  solved  to  yield  the  distribution 
function  of  electrons  in  action  space.  Normalizing  the 
distribution  function  to  the  globally  stochastic  fraction 
of  the  equilibrium  electrons  provides  a  leakage  current. 
The  analysis  yields  veiy  good  agreement  with  the  exper¬ 
imental  results  of  Orzechowski  and  Bekefi  (2)  on  leakage 
currents  in  a  magnetically  insulated  device. 


Equilibria 

The  equilibria  used  in  the  analysis  consist  of  a  constant- 
density,  non-relativistic  electron  sheath,  a  uniform  ap¬ 
plied  magnetic  field  B  =  B0i,  and  an  applied  voltage 
V  across  a  planar,  periodic  gap  of  length  L  and  width 
d.  The  equilibrium  sheath  is  assumed  to  be  confined  on 
the  cathode  side  with  height  x.,  where  0  <  x.  <  d.  The 
geometry  is  illustrated  in  Fig.  1.  The  i  direction  is  ig¬ 
nored  throughout  the  analysis.  The  mean  flow  is  in  the 
y  direction. 


Anode 

®8o 


Sheath 


Cathode 


T 

-T  0 
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Fig.  1.  Geometry  and  coordinate  system. 


The  single  particle  Hamiltonian  for  the  region  0  < 
x  <  x.  can  be  written  in  the  dimensionless  form 


"  2  +  2  **  2  ’ 


C) 


where  Q  =  (2e^./mfl*xJ)  and  4>.  is  the  potential  at 
the  sheath  surface.  Distances  are  normalised  to  x.,  mo¬ 
menta  to  mflx.,  and  time  to  1/(7.  The  parameter  Q 


is  equal  to  ui2jU2,  where  uip  is  the  electron  plasma  fre¬ 
quency  and  fl  is  the  electron  cyclotron  frequency.  The 
Hamiltonian  (I)  defines  a  class  of  systems  parametrized 
by  Q.  The  well-known  Orillouin  flow  equilibrium  (3j  is 
obtained  with  Q  —  I.  Equilibria  exist  for  the  range 


0  <  Q  <  1. 


The  equilibrium  Hamiltonian  can  be  cast  in  conve¬ 
nient  form  for  perturbation  analysis  through  the  use  of 
canonical  transformations.  The  Hamiltonian  in  action- 
angle  variables  |l|  becomes 


where 


ui'2  =  1  -  Q,  and  J+  =  (L/2n)Pt.  The  two  character¬ 
istic  frequencies  are  ug  —  w,  which  corresponds  to  the 
^-direction  bounce  frequency,  and  tjg  =  (2x/L)JJ^(l  - 
which  corresponds  to  the  y-direction  transit  fre- 
qncnc.y. 


Stochasticlty  and  Diffusion 
In  the  absence  of  perturbations,  an  electron  initially  con¬ 
fined  to  the  sheath  region  may  remain  so  indefinitely. 
This  permits  complete  insulation  for  a  sufficient  mag¬ 
netic  field.  For  non- relativistic  systems  in  which  all  the 
electrons  have  zero  energy  //  and  canonical  momentum 
P„  this  critical  field  is  given  by  £JC  =  c(2mV/e)*/J/d 
[4|.  In  the  presence  of  perturbations,  the  situation  is 
much  more  complicated.  For  small  perturbations,  local 
invariants,  known  as  KAM  surfaces  |5-7|,  may  persist 
that  serve  to  barricade  the  electron  motion  and  prevent 
leakage.  The  analysis  in  our  previous  paper  was  devoted 
primarily  to  establishing  the  conditions  under  which  the 
KAM  surfaces  are  destroyed  and  electron  leakage  be¬ 
comes  possible.  We  considered  a  dimensionless  pertur¬ 
bation  potential  of  the  form 

^  =  A  ^  z(  1  -  x)  sin(2 wy/L),  (4) 

where  A  is  an  arbitrary  amplitude  factor  and  L  m  the 
wavelength  of  the  perturbation  normalised  to  The 


perturbation  electric  field  and  the  background  electric 
field  are  related  in  magnitude  by  \EVm^JEx.\  =  2A/Z, 
The  use  of  action-angle  variables  allows  us  to  write  the 
perturbed  Hamiltonian  in  the  convenient  form 


OO 

tf  =  ff0  +  A  £  Hk(J)ain(<f>-kO),  (5) 

Jfc=  — OO 

where  Ho  is  given  by  Eq.  (2).  We  refer  the  reader  to 
Ref.  1  for  details  on  the  functional  form  of  the  Hk.  The 
Hamiltonian  (5)  exhibits  resonances  at  the  stationary 
phase  points  of  the  driving  term,  i.e.,  where  j>  -  k$  =  0. 
The  frequencies  are  approximately  given  by  the  equi¬ 
librium  Hamiltonian.  In  regions  of  phase  space  where 
the  equilibrium  frequencies  satisfy  the  resonance  condi¬ 
tion,  resonant  islands  will  appear.  Roughly  speaking,  if 
the  islands  are  large  enough  to  overlap,  the  motion  of 
particles  near  the  islands  becomes  very  irregular.  This 
is  known  as  intrinsic  stochasticity.  The  analysis  in  our 
previous  paper  indicated  that  the  stochasticity  is  a  sensi¬ 
tive  function  of  Q.  The  level  of  stochasticity  is  depicted 
by  the  stochasticity  parameter  S,  which  represents  the 
fraction  of  intervening  phase  space  covered  by  adjacent 
resonances.  An  approximate  criterion  for  stochasticity 
is  due  to  Chirikov  [8|  and  is  given  by  S  >  I.  A  rough 
approximation  to  the  stochasticity  parameter  calculated 
in  our  previous  paper  is 


S  < 


8  Q  (*&\>/, 


(«) 


The  largest  values  of  5  are  found  on  the  II  =  0  energy 
surface,  where  the  relation  approaches  equality.  Note 
that  S  diverges  as  Q  — ►  1:  the  Brillouin  flow  limit.  The 
Brilkmin  flow  equilibrium  is  globally  stochastic  for  arbi¬ 
trarily  small  perturbations. 

In  Ref.  I  we  arrived  at  the  following  qualitative  pic¬ 
ture  of  diode  saturation.  The  electron  sheath  diffuses 
towards  the  anode  as  a  result  of  stochastic  motion.  The 
parameter  Q  decreases  as  the  density  profile  spreads  out 
over  the  gap.  If  the  perturbation  is  sufficiently  large, 
the  sheath  will  spread  to  the  anode  (r,  — >  d)  and  a  leak¬ 
age  current  will  be  established.  The  assumption  that 


the  density  profile  remains  relatively  flat  yields  a  sim¬ 
ple  relation  between  the  saturated  value  of  Q  and  flc, 
namely,  Q,K,  —  (FI,/ FI)1  Hence,  the  larger  the  value  of 
/l//Jr,  the  more  difficult  electron  leakage  becomes.  In 
the  following  analysis,  we  proceed  to  quantify  this  leak¬ 
age  current. 

A  quasilinear  diffusion  coefficient  for  the  stochasiti- 
cally  wandering  electrons  is  obtained  with  the  standard 
technique  of  integrating  the  perturbation  force  over  the 
unperturbed  trajectories  to  calculate  the  change  in  ac¬ 
tion  for  a  given  initial  condition  as  a  function  of  time. 
Averaging  over  the  initial  phases  allows  us  to  form 
D(.Ig;H).  In  the  interest  of  brevity,  we  omit  the  anal¬ 
ysis  here;  details  may  be  found  in  Ref.  10.  The  diffu¬ 
sion  coefficient  is  given  by  D(J§\H)  =  7s(Ak//,)1/4, 
where  7#  =  2n/u>g  and  *(./»)  =  u+(Jg)/ut.  The  dif¬ 
fusion  coefficient  D(Jg\  If)  is  found  to  be  a  monoton- 
irally  increasing  function  of  J$.  This  implies  a  mean 
shift  C  =  dD/dJg  [9|  as  well  as  a  spreading  of  the  single 
particle  distribution  function.  Since  the  perturbation  we 
consider  is  time-independent,  the  motion  of  a  given  elec¬ 
tron  is  confined  to  a  constant  energy  surface.  Therefore, 
each  energy  surface  has  its  own  diffusion  equation. 

Steady-State  Diffusion 

Having  obtained  the  diffusion  coefficient,  we  consider  the 
diffusion  equation  for  the  stochastic  electrons: 

where  S(./»;  H)  is  a  source  term.  In  our  system,  the  only 
source  of  electrons  is  the  cathode  at  (Jg,  H)  =  (0, 0). 
Since  the  stochastic  motion  is  confined  to  constant  en¬ 
ergy  surfaces,  the  steady-state  solution  will  be  non-zero 
only  along  FI  =  0.  The  other  stochastic  surfaces,  which 
are  without  a  source,  will  become  depleted  of  electrons. 
Along  II  =  0,  we  can  replace  the  source  by  rS(Jg). 
Setting  df/dt  =  0  and  integrating  over  J$,  we  obtain 
the  well-known  flux  relation  T  =  -D(J$)df /9J$  - 
constant.  We  may  integrate  once  again  to  obtain  the  dis¬ 
tribution  function.  The  anode  is  modeled  as  a  perfect 


absorber  of  electrons;  the  boundary  condition  is  then 
/(Ja)  =  0,  where  Ja  is  the  action  value  correspond¬ 
ing  to  the  anode.  Graphing  D(J$)  on  a  log-log  plot 
indicates  that  D(J$)  may  be  accurately  fitted  by  a  sim¬ 
ple  power-law  expression:  D(J»)  =  Ai~i(Q,  L)Jgv^'L\ 
where  q(Q,  L)  and  i>(Q,  L)  contain  the  parametrization 
of  D(Jg)  in  terms  of  the  equilibrium  and  wavelength 
parameters.  For  a  wide  range  of  Q  and  L  we  find 
1.6  <  u  ;$  1.8.  The  integration  of  the  flux  equation 
is  now  straightforward: 

nJ,)  -  j 11  -  "I  •  '•» 

Note  that  the  distribution  is  proportional  to  the  fiux  T, 
which  is  yet  to  be  determined. 

Leakage  Currents 

Integrating  Eq.  (8)  from  J#  =  0  to  J§  -  J a  relates 
the  flux  T  to  the  number  of  diffusing  electrons:  T  = 
(2  -  v)D(J a)J a~* Pi  where  p  is  defined  by  the  integral 
of  f(J$)  from  0  to  Ja-  To  calculate  the  actual  leakage 
current,  we  must  determine  p.  The  diffusing  population 
p  is  required  to  be  equal  to  the  globally  stochastic  frac¬ 
tion  of  the  equilibrium  electrons. 

The  calculation  of  p  consists  of  a  numerical  inte¬ 
gration  over  the  allowed  energy  surfaces.  A  contribu¬ 
tion  equsd  to  the  local  equilibrium  density  is  made  if  a 
given  energy  surface  is  globally  stochastic.  The  maxi¬ 
mum  value  of  p  in  our  dimensionless  units  is  1.0,  corre¬ 
sponding  to  global  stochasticity  on  all  available  energy 
surfaces.  Details  on  the  calculation  of  p  can  be  found 
in  Ref.  10.  For  comparison  with  earlier  experimental  re¬ 
sults  and  to  obtain  a  useful  measure  of  the  leakage  cur¬ 
rent,  we  normalise  the  leakage  current  er  to  the  Child- 
Langmuir  |1 1]  current  for  the  actual  gap  d.  We  obtain 

£  =  (9) 

The  plot  in  Fig.  2  compares  our  theoretical  results  (solid 
curve)  to  a  simple  fit  to  the  experimental  data  of  One- 
chowski  and  Bekefi  (dashed  line).  The  circles  are  the 


Fig.  2.  Leakage  currents  I/lti  versus  D/Bc 
f«*r  A  —  0.065  ami  L  —  2.7.  The  dashed  line  is  a 
fit.  to  the  experimental  results  of  Ref.  2  for  their 
mechanical  gap  0.5-1  cm.  The  dotted  line  is  the 
theoretical  prediction  of  Ref.  12.  The  circles 
are  the  result  of  a  semi-numerical  calculation. 

result,  of  a  seini-numerical  calculation  in  which  we  deter¬ 
mine  the  mean  inverse  crossing  time  (rj 1 )  numerically 
to  yield  the  semi-numerical  flux  rBllnl  =  p(r^ 1 ).  Also 
shown  is  the  result  of  Mouthaan  and  Susskind  (12|  (dot¬ 
ted  line),  an  earlier  attempt  to  calculate  leakage  cur¬ 
rents  in  r.rosscd-field  devices.  The  agreement  between 
both  our  analytical  and  semi-numerical  calculations  and 
the  experimental  results  is  very  good.  Similar  curves  are 
obtained  that  fit  the  remaining  data  of  Oreechowski  and 
Hekefi.  Comparison  of  our  results  with  their  experimen¬ 
tal  data  indicates  that  the  characteristic  wavelength  of 
perturbations,  normalised  to  the  actual  gap  d,  was  an 
increasing  function  of  the  mechanical  gap.  Had  the  per¬ 
turbations  been  a  result  of  mechanical  irregularities,  the 
converse  would  have  been  true.  For  their  mechanical  gap 
of  0.54  cm,  the  characteristic  wavelength,  normalised 
to  the  actual  (dynamical)  gap,  was  approximately  2.7; 
the  electric  Held  perturbation  was  approximately  4.8%  of 
the  accelerating  field.  For  their  mechanical  gap  of  0.23 
cm,  the  normalised  wavelength  was  approximately  2.2; 
the  corresponding  transverse  electric  field  was  approxi¬ 
mately  7.3%  of  the  accelerating  Held. 

It  is  clear  from  examining  the  theoretical  curve  that 
this  analysis  predicts  a  sharp  cutoff  at  some  value  of 
D/Dt  greater  than  I  (1.3  in  Fig.  2).  The  location  of  the 


cutoff  can  be  calculated  using  the  approximate  relation 
for  the  stochasticity  parameter  given  by  Eq.  (6).  Setting 
S  equal  to  1  and  replacing  Q  with  (B,c/Bc)7,  we  find 


This  relation  defines  a  revised  critical-insulating-field  for 
magnetic  insulation,  which  we  call  “B-supcrcritical"  and 
denote  by  BK. 

It  is  important  to  note  that  the  above  analysis  per¬ 
tains  to  the  actual  or  dynamical  gap  d,  whereas  the  au¬ 
thors  of  Ref.  2  calculated  B/Bc  with  the  claim  that  the 
actual  gap  was  nearly  the  mechanical  gap.  Recent  ex¬ 
perimental  work  [I3|  has  initiated  a  quantitative  inves¬ 
tigation  of  the  actual  gap;  their  results  indicate  that 
the  expansion  of  the  electrode  plasmas  may  significantly 
reduce  the  actual  gap.  With  an  understanding  of  the 
actual-gap  dynamics,  more  accurate  tests  of  this  theory 
may  be  attempted. 


Acknow  ledgments 

This  work  was  supported  by  DOE  Contract  No.  DE- 
AS08-81DP40139  and  by  ONR  Contract  No.  N00014- 
85-K0212. 


References 

1.  M.  P.  Desjarlais  and  R.  N.  Sudan,  Phys.  Fluids  20, 
1245  (1986). 

2.  T.  J.  Orxechowski  and  C.  Bekefi,  Phys.  Fluids  10, 
43  (1976);  23,  978  (1979). 

3.  L.  Brillouin,  Phys.  Rev.  67,  260  (1945). 

4.  A.  W.  Hull,  Phys.  Rev.  18,  31  (1921). 

5.  A.  N.  Kolmogorov,  Akad.  Nayk  SSSR  Doklady  08, 
527  ( 1954)  (in  Stochastic  Behavior  in  Classical  and 
Quantvm  Hamiltonian  Systems,  edited  by  G.  Casati 
and  J.  Ford  (Springer,  New  York,  1979),  p.  5 1). 

6.  V.  I.  Arnold,  Russ.  Math.  Surv.  18,  9,  85  (1963). 

7.  J.  Moser,  Nachr.  Akad.  Wise.  Gottingen  Math. 
Phys.  Kl.  3,  1  (1962). 

8.  B.  V.  Chirikov,  Phys.  Reports  52,  265  (1979). 

9.  L.  D.  Landau,  Zh.  Eksper.  Theor.  Fis.  7, 203  (1937). 

10.  M.  P.  Desjarlais,  Ph.  D.  thesis,  Cornell  University 
(to  be  submitted). 

11.  C.  D.  Child,  Phys.  Rev.  32,  492  (191 1);  I.  Langmuir, 
Phys.  Rev.  3,  450  (1913). 

12.  K.  Mouthaan  and  C.  Susskind,  J.  Appl.  Phys.  37, 
2598  (1966). 

13.  Y.  Maron,  M.  D.  Coleman,  D.  A.  Hammer,  and  H.- 
S.  Peng,  these  proceedings. 


-5- 

ION  BEAM  OPENING  SWITCH* 

J.B.  Greenly,  g.D.  Rondeau, 

H.T.  Sheldon,  D.A.  Hanmer,  P.L.  Dreike** 

Laboratory  of  Plasma  studies 
Cornell  University 
Ithaca,  N.y.  14853 


Abstract 

Plasma  opening  switches  (POS)  have  shown 
excellent  characteristics  in  pulsed-power 
applications.  Proposed  POS  scaling  predicts 
that  the  fastest  opening  time  for  a  given 
conducted  current  should  occur  using  a  high- 
velocity,  low  density  plasma  as  the  switch  med¬ 
ium.  The  Ion  beam  opening  switch  (IBOS)  uses  a 
charge-neutral  ion  beam  of  100-300  kv,  <120  A/cm2 
as  the  switch  "plasma".  Its  velocity  of  200-600 
cm/us  and  density  of  -1012/cm3  make  this  a  very 
fast,  low-density  plasma  compared  with  typical 
-10  cm/usec  and  ~1013/cm3  POS  plasmas.  The 
IBOS  has  conducted  >70  kA  flowing  in  a  parallel 
plate  transmission  line  driven  by  a  4  fl  pulser. 
IBOS  opening  time  is  load  dependent  being  4  4  ns 
into  a  15  nH  and  about  twice  as  long  into  an 
electron  diode  load.  However,  switch  inpedance 
is  not  zero  while  current  is  carried,  and  rise  to 
<3o  near  peak  current.  Current  conducted  before 
opening  is  not  linear  with  either  injected  ion 
current  or  switch  area.  The  results  suggest  that 
electron  current  conduction  occurs  in  a  narrow 
channel  of  $1  cm  axial  length.  In  this  experi¬ 
mental  geometry,  replacing  the  IBOS  with  a 
standard  POS  gave  essetentially  identical  be¬ 
havior,  except:  first,  an  order  of  magnitude  less 
injected  flux  was  required  to  conduct  a  given 
current  with  the  POS  than  with  IBOS,  correspond¬ 
ing  to  similar  injected  plasma  densities  for  the 
two  types  of  switches.  Second,  clear  evidence  of 
i  *  p  motion  of  the  POS  current  was  seen,  while 
none  was  evident  with  IBOS.  We  are  developing  a 
model  of  this  switch  )r;  which  electrons  conduct 
current  by  6  X  0  drift,  that  is  a  9elf-consistent 
Hall  current  in  the  self-magnetic  field  of  the 
line. 

Introduction 

Plaana  opening  switches  (POS)  have  been 
shewn  to  be  of  great  interest  to  the  continued 
development  of  pulsed  power (1| .  In  the  shortest 
pulse  application  these  switches  are  Intended  to 
carry  large  currents  (>1  MA)  and  then  open  in 
times  of  order  10ns  or  less.  They  may  be  able  to 


provide  large  gains  in  voltage  and  power  into  a 
load  ty  opening  while  current  is  flowing  in  an 
inductance  "charged"  by  a  standard  pulsed-power 
generator.  A  POS  is  made  by  injecting  a  flowing 
plasma  through  one  (perforated)  electrode  usually 
the  anode  of  a  magnetically-insulated  vacuun 
transmission  line,  usually  in  coaxial  geometry. 
When  a  high  power  pulse  propagating  down  the  line 
arrives  at  the  location  of  the  plasma  the  current 
is  initially  short  circuited  throucfi  the  plasma 
until  the  current  reaches  a  certain  level,  at 
which  time  the  impedance  of  the  plasma  suddenly 
rises  to  a  high  value  (i.e.  the  switch  "opens"). 

A  model  based  on  plasma  erosion,  or  the  sweeping 
out  of  plasma  from  the  gap  of  the  transnission 
line  as  the  current  increases,  has  been 
developed.  This  model,  called  the  plasma  erosion 
opening  switch  (PEGS) 111,  predicts  opening  of  the 
switch  to  high  impedance  as  magnetic  insulation 
is  established.  This  model  predicts  that  the 
current  carried  by  the  switch  plasma,  before 
opening  begins  and  voltage  appears  across  the 
line,  is  proportional  to  the  flux  of  plasma  ions 
reaching  the  cathode.  Also,  the  opening  rate, 
once  this  "conduction"  current  is 
exceeded,  is  inversely  related  to  the  plasma 
density.  Thus,  to  make  an  effective  POS,  the 
PB0S  model  determines  the  required  plasma  flux 
for  a  given  current  to  be  conducted.  Further¬ 
more,  it  indicates  that  fastest  opening  should 
occur  when  that  flux  is  provided  by  a  high 
velocity,  lew-density  injected  plasma. 

Previous  POS  experiments  have  typically  used 
plaana  sources  of  -1013  on3  density  and 
10-20  an/us  velocity.  In  order  to  test  the  PEGS 
scaling  predictions,  we  have  performed  an 
experiment  in  which  the  injected  plasma  is  a 
nearly  charge  -  and  current-neutral  ion  beam, 
usually  protons  of  >  400  arvAis  velocity  (>100 
keV)  and  about  *10l2/cm3  density. 

A  second  reason  for  interest  in  this 
experiment  is  the  expectation  that  such  high 
plaana  ion  velocities  should  eliminate  J  x  g 
motion  of  the  witch  plasma  during  conduction 
since  the  ion  larmor  radius  is  always  larger  than 
the  line  gap. 


Description  of  the  Experiment 

The  experiment,  called  an  IBOS,  for  ion  beam 
opening  switch,  is  shewn  in  Fig.  1.  A  flat 
parallel-plate  transmission  line,  10  cm  wide  with 
a  1  cm  gap,  was  pulsed  in  negative  polarity  so 
the  ion  beam  from  a  simple  planar  magnetically 
insulated  ion  diode  (3]  could  be  injected  through 
a  10  cm  x  10  cm  area  of  the  grounded  line  anode. 
The  line  was  magnetically  insulated  at  <  10  kA, 
to  the  highest  voltage,  -450  kv,  applied.  The 
inductance  between  the  output  of  the  4n  ONI  II 
water  pulse-forming  line  and  the  switch  region 
was  -ISOnH,  while  the  load  inductance  was  *15  nH. 
Up  to  12  kA  of  injected  ion  flux  was  measured  by 
the  array  of  15  biased  charged  collectors  mounted 
in  the  switch  cathode. 

The  bogm  uniformity  was  better  than  t  30%  over 
the  switch  voa.  The  ion  bean  pulse  lasted  >200 
ns,  and  the  line  current  arrival  time  oould  be 
varied  throughout  this  pulse,  lb  compare  the 
properties  of  IBOS  with  a  typical  POS,  the  ion 
diode  was  replaced  by  a  standard  Mendel-type  (2) 
plasma  gun  placed  20  cm  from  the  line  cathode. 

Experiment  Results 

The  ion  beam  did  produce  an  excellent 
opening  switch  into  diode  loads.  Figure  2  shows 
an  IBOS  shot  using  a  blunt  carbon-cathode 
electron  diode  load  of  about  4n.  the  load 
inpedance  was  infinite  (electron  emission  did  not 
begin)  until  the  switch,  and  load,  voltage 
reached  *150  kV.  The  downstream  (load)  current 
then  rose  rapidly  in  about  -10  n3,  while  the  up¬ 
stream  current  in  the  line  inductance  dropped 
somewhat,  adding  an  LI  voltage  of  <.100  kV  at 
the  switch.  Almost  all  of  the  current  was 
diverted  to  the  load  after  opening.  It  is 
notable  that  the  anount  of  current  conducted  with 
negligable  voltage  across  the  switch  was  quite 
small,  perhaps  10”15  kA,  although  in  this  shot 
the  injected  ion  flux  was  11-12  kA.  The  PE0S 
model  would  predict  a  conduction  current  of  at 
least  40  times  the  ion  flux  (if  the  ions  are  all 
protons).  This  result  was  universally  observed 
with  I BOR  shots:  the  conduction  current  before 
switch  voltage  appeared  in  excess  of  the  -10  kV 

inductive  voltage  due  to  I  in  the  switch  was 
quite  <snall,  and  far  below  tlie  bipolar-flow  ratio 
characteristic  of  PEGS  model  behavior  and  pre¬ 
vious  PCS  experimental  results.  When  IBOS  was 


used  with  a  short-circuit  ( 15nH)  load,  load  cur¬ 
rent  began  very  early  and  very  little  switching 
could  be  produced,  as  would  be  expected  since 
load  voltage  appears  so  early.  These  results 
were  anticipated,  since  bipolar  flow  with  200kV 
injected  ions  would  require  200  kv  counter  flow¬ 
ing  electrons  in  the  switch  gap;  thus  P0O5  type 
"conduction"  current  would  not  be  reached  until 
the  switch  voltage  approached  200  kV.  Even  at 
this  voltage  however,  the  observed  switch  current 
was  well  below  the  bipolar  value  except  when  the 
switch  axial  length  was  reduced  to  <  1  cm,  as 
shewn  belcw. 

Wien  we  replaced  the  ion  beam  with  the  plasma 
gun,  our  POS  results  did  agree  qualitatively  with 
previous  experiments  in  terms  of  ratio  of  (low 
voltage)  conduction  current  to  plasma  ion  flux. 
Figure  3  shows  a  POS  shot  into  the  short-circuit 
( 15  nH)  load  in  which  a  conduction  current  of 
about  25  kA  was  carried  with  a  plasma  flux  of 
about  0.3  kA. 

Since  IBOS  nevertheless  gave  excellent 
switching  into  a  typical  diode  load  which  did  not 
"close"  and  begin  to  carry  current  until  a 
certain  voltage  was  reached,  we  investigated  this 
behavior  further.  In  order  to  simplify  the 
effect  of  the  load  on  the  switch-circuit 
behavior,  we  made  a  series  of  shots  in  which  the 
electron  diode  was  replaced  by  the  15  nH 
short-circuit  with  an  added  series  surface- 
flashover  switch  which  broke  down  at  about  200 
kV.  Thus,  the  load  began  as  an  open  circuit,  as 
does  a  diode  load,  but  after  this  load  "turned 
on",  it  had  a  simple  15  nH  impedance  instead  of 
the  more  carpi icated  time-decreasing  inpedance  of 
a  diode  load.  A  typical  IBOS  shot  with  this  load 
is  show:  in  Fig.  4.  The  switch  (upstream) 
current  and  voltage  rise  as  with  the  diode  load 
until  the  surface  switch  closes.  Then  the  load 

current  rises  very  rapidly,  with  full  current 
diverted  to  the  load  in  as  little  as  4  ns. 

Again,  a  negative  I  can  appear  upstream  during 
opening,  generating  inductive  voltage  which 
appears  as  the  plateau  in  Vg*.  After  opening, 
the  voltage  drops  to  the  low  level  characteristic 
of  L  I  into  the  15  rH  load. 

When  a  POS  was  used  with  the  sane  load. 

Fig.  5,  the  switching  waveforms  were 
indistinguishable  from  IBOS.  However,  the  plasma 
flux  required  for  a  given  current  conducted 
before  the  surface-switch  breakdown  voltage  was 
reached  was  again  an  order  of  magnitude  less  for 


PQS  than  for  IBOS.  With  this  load  the  scaling  of 
current  at  the  time  of  opening  with  IBOS  and  POS 
flux  was  investigated.  Results  are  shown  in 
Fig.  6.  The  switch  current  did  not  scale 
linearly  with  flux  for  either  tbos  or  POS.  Both 
scale  more  nearly  with  the  square  root  of  the 
flux.  Although  the  flux  differed  greatly  between 
IBOS  and  POS  for  the  same  switch  behavior,  the 
density  as  determined  from  measured  flux  and 
velocity  for  similar  IBOS  and  POS  shots,  was 
often  within  a  factor  of  two  of  agreement. 

HP  also  investigated  the  scaling  of  IBOS  and 
POS  switch  current  with  the  axial  length  of  the 
switch,  by  maskinq  off  various  parts  of  the  anode 
extrance  area,  keeping  the  incident  plasma  or 
beam  flux  density  constant.  Results  are  shewn  in 
Fig.  7.  In  both  IBOS  and  POS,  the  current 
carried  before  flashover  switch  voltage  was 
reached  was  far  from  proportional  to  switch 
length  (area)  except  for  very  short  switches. 

One  possible  interpretation  of  this  result  is 
that  the  switch  current  is  being  carried  in  an 
axially  very  narrow  channel,  so  when  the  switch 
plasma  axial  length  is  greater  than  this  channel 
thickness,  very  little  additional  current  is 
carried  by  the  rest  of  the  plasma.  The  experi¬ 
mental  results  would  inply  current  channel  axial 
length  of  the  order  of  1/2-lcm.  If  this  is  true, 
a  two  dimensional  model  of  the  switch  behavior  is 
certainly  required,  and  the  simplest  PEGS  model 
cannot,  entirely  b°  expected  to  explain  the 
behavior  of  these  switches. 

Wa  are  developing  a  new  opening  switcli  model 
based  on  a  self-consistent  §  x  9  Hall  current 
channel  cross  the  switch  plasma  gap.  This 
channel  results  from  the  electric  field  of  an 
axial  charge  separation  due  to  V  x  B  force  on  the 
current  carrying  electrons,  and  the  self-magnetic 
field  of  that  current.  A  cathode  sheath  is 
required  to  Inject  electrons  at  the  self  con¬ 
sistent  drift  velocity  Into  the  Hall  channel, 
giving  a  relation  between  the  switch  (sheath) 
voltage  and  switch  current.  In  particular,  if 
the  Hall  channel  thickness  is  c/wpe,  the 
collisionless  skin  depth,  the  switch  (sheath) 
voltage  becomes  equal  to  the  axial  Hall  voltage, 
which  varies  inversely  with  the  plasma  density 
and  proportional  to  switch  current  squared.  Thus 
at  a  given  switch  voltage,  e.g.  the  voltage  at 
which  the  surface  switch  broke  down  In  our 
experiment,  this  model  predicts  switch  current 
proportional  to  the  square  root  of  the  pla.ma 


density.  The  collisionless  skin  depth  in 
experiments  was  <  0.5  cm.  A  detailed  presen¬ 
tation  of  this  model  will  appear  in  a  future 
publication.  This  model  gives  switch  iirpedance 
history,  and  scaling  with  switch  length  and 
plasma  parameters,  that  agree  with  both  our  IBOS 
and  POS  experiments. 

One  further  difference  between  IBOS  and  POS 
results  was  observed  which  agreed  with  the 
expectation  that  IBOS  might  be  iirmune  to  J  x  p 
motion.  No  damage  was  ever  observed  beyond  the 
downstream  boundary  of  the  beam  entrance  area 
with  IBOS,  while  with  the  POS,  one  shot  sufficed 
to  melt  a  mylar  sheet  defining  the  anode  entrance 
aperture  for  a  distance  of  2  -  4  cm  past  the 
downstream  edge  of  the  switch  region.  This 
damage  is  presumably  due  to  current-carrying 
electrons  accelerated  across  the  switch  before 
and  during  opening.  An  estimate  of  the  distance 
the  POS  plasma  should  be  pushed  downstream  by 
J  *  6  indicates  that  if  the  current  were  being 
carried  by  the  full  mas3  of  the  switch  plasma, 
very  little  motion  should  occur,  whereas  if  the 
current  channel  were  only  1  cm  long,  motion 
should  be  comparable  to  that  observed.  However, 
sane  "pinching"  of  the  current  in  the  direction 
across  the  10  cm  width  of  the  line  is  also 
evident,  which  might  complicate  an  estimation  of 
the  expected  motion. 
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Fig.  1  (A)  OMNI  It,  4a,  kV,  100  ns;  interface  to  (B) 
vacuum  chamber,  driving  (C)  switch  line  cathode  and 
(D)  anode,  10  cm  wide  plates.  Diagnostics:  (E) 
upstream  and  (F)  downstream  current;  (G)  inductive 
divider  for  switch  voltage;  (H)  multi-aperture  ion 
beam  current  monitor.  Load  (I)  surface-f lashover 
series  switch;  (J)  electron  diode.  Ion  diode: 
cathode  (K) ;  anode  (L)  driven  SNLA  Marx,  300  kV,  10n, 
magnetically  insulated  by  (M)  capacitor  bank. 
Flux-excluding  tube  (N)  tr-<•’=rv"'^,■<,  ion  beam  (P). 


225kV 

35  k  A 

35kA 

20  ns  /  div 

Fig.  2.  IBOS  switching  into  an 
electron  diode  load.  Vsw  is  the 
inductive  voltage  divider,  Iu  is 
sun  of  upstream  B  loops,  and  I<j 
is  sum  of  downstream  B  loops. 


20  ns /div  20  ns /div  20ns/div 


Fid.  3.  POS  switching  into  a 
15  nfl  load. 


Fig.  4  190S  switching  into  15  nil 
plus  series  surface  swU-ch. 


Fig.  5  PtJS  switching  into  15  nH 
plus  series  surface  switch. 


Fig.  6  Switch  current  at  opening  versus  injected 
flux  for  POS  and  IBPR. 


SwMtfl  t» 


Fig.  7  Switch  current  (Igw*  Iu)  at  opening  versus 
switch  axial  length  for  POS  and  IBOS. 
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Abstract 

An  active  anode  plasma  source  has  been 
operated  on  the  LONGSHOT  annular  magnetically 
insulated  ion  diode.  This  source  uses  an 
inductive  voltage  from  a  single  turn  coil  to 
break  down  an  annular  gas  puff  produced  by  a 
supersonic  nozzle.  The  resulting  plasma  is 
magnetically  driven  toward  the  radial  insulating 
magnetic  field  in  the  diode  accelerating  gap  and 
stagnates  at  a  well-defined  surface  after  about 
300  ns  to  form  a  plasma  anode  layer  defined  by 
magnetic  fields.  An  ion  beam  is  then  extracted 
from  this  plasma  layer  ty  applying  typically  a 
150kV,  1  us  pulse  to  the  accelerating  gap. 
Optimization  of  the  timing  of  the  gas  puff,  the 
plasma  production  discharge  and  the  high  voltage 
pulse  has  resulted  in  1  us  duration  75-150  keV 
ion  beam  pulses  with  >100  A/cm2  peak  ion  current 
density  over  an  area  of  about  400  an2.  Up  to 
5J/cm2  has  been  collected  by  a  4  cm2  callmeter. 
The  diode  inpedance  history  can  be  varied  so  that 
rising,  flat,  and  falling  voltage  pulse  waveforms 
can  be  produced.  Streak  photographs  of  beamlets 
inpinging  on  a  scintillator  and  time  integrated 
targets  both  show  beam  divergence  angles  <  3°, 
but  under  certain  operating  conditions,  large 
excursions  (-25°)  in  mean  aiming  angle  on  time 
scales  of  20-200  ns.  These  and  other  operating 
characteristics  of  the  gas-beakdcwn  diode  are 
discussed. 

Introduction 

Magnetically  insulated  diodes  (MID's)  have 
proven  to  be  very  useful  intense  ion  beam  sources 
in  applications  ranging  frcm  Inertial  Confinement 
Fusion  at  the  >1012W  level  (1)  to  materials 
science  research  at  a  power  level  of  about  10loW 
[2| .  The  ion  source  at  the  anode  in  these  diodes 
has  usually  been  a  flashover  plasma  induced  on  a 
dielectric  surface  ty  the  high  voltage  pulse  that 
accelerates  the  ion  bean  [3] .  Electrons  which 
"leak"  across  the  insulating  magnetic  field  and 
inplnge  upon  the  surface  are  also  believed  to 
play  a  major  role  (4).  Such  "surface  flashover 


anodes"  have  several  disadvantages,  including: 

1.  The  turn-on  delay,  especially  for  short 
pulses  or  at  low  voltages,  can  waste  a  sub¬ 
stantial  fraction  of  the  power  pulse  delivered  to 
the  ion  diode?  2.  The  large  nunber  of  neutrals 
blown  off  the  surface  as  part  of  the  flashover 
process  are  believed  to  be  responsible  for  the 
rapidly  inpedance  characteristic  of  these  diodes 
(7)?  3.  Ion  beams  drawn  from  surface  flashover 
plasmas  tend  to  be  a  mixture  of  protons  and 
heavier  species,  especially  various  charge  states 
of  carbon  and  oxygen  and  4.  The  life  of  surface 
flashover  anodes  ranges  frcm  one  pulse  at  the 
highest  power  levels  [11  to  perhaps  a  few  hundred 
at  the  lO’w/cm2  level  (2,7)  due  to  the  damage 
they  sustain  from  electron  borrtxardment. 

There  have  been  several  attempts  to  provide 
an  anode  plasna  ion  source  which  does  not  suffer 
frcm  these  disadvantages,  such  as  the  plasma 
filled  diode  of  Mendel  [8],  the  plasma  soiree 
used  by  Hunphries  et  al  [9)  which  is  a  direct 
ancestor  of  the  source  we  describe  here,  and  the 
actively  driven  surface  flashover  anodes  (using 
ari  energy  source  independent  of  the  main  ion 
accelerating  power  pulse)  of  McClure  [101  and 
Greenly  [111.  Other  possibilities  in  various 
stages  of  development  are  described  in  other 
papers  in  this  conference  [O-E-6,  O-H-2  and  4, 
P-B-20,  etc) . 

In  this  paper  we  describe  a  magnetically 
insulated  diode  with  a  gas-breakdown  anode  plasna 
ion  source  which  improves  upon  almost  all 
categories  of  results  obtained  with  a  surface 
flashover  anode  when  operated  on  the  same  pulsed 
power  generator  (LONGSHOT).  The  total  icn  output 
has  been  more  than  doubled.  A  variety  of  diode 
inpedance  histories  including  constant  or  even 
rising  inpedance  can  be  selected  by  changing 
parameters  of  the  source.  Proton  beams  produced 
when  Hj  gas  is  used  show  no  inpurity  species, 
within  the  t  10%  ability  of  the  foiled  Faraday- 
cup  measurement  technique.  Furthermore,  the  beam 
extraction  can  be  made  to  begin  coincident  with 
the  beginning  of  the  voltage  pulse.  The  diode 
was  fired  hundreds  of  times  without  damage  or 


replacement  of  parts.  Oily  the  divergence  angle, 
3°,  was  not  improved  over  surface  flashover  diode 
results. 

Clear  correlation  is  seen  between  various 
aspects  of  the  diode  performance  and  measured 
anode  plasma  ch->>  ’-'teristics.  For  example,  the 
ion  beam  current  density  is  related  to  the  anode 
plasma  flux  and  to  its  position  as  determined  by 
the  magnetic  field  of  the  plaiana  source  and 
diode.  These  dependences,  along  with  the  ion 
beam  chacteristics,  have  led  to  an  understanding 
of  the  source  gap  operation  which  will  be 
dismissed  in  future  pel  1  ications. 

Description  of  the  Gas-Breakdown  Anode  Plasma 
Source. 

The  principal  conponents  of  the  gas-break¬ 
down  anode  plasma  source  are  shown  in  the  LONG¬ 
S’ lOT  diode  in  Fig.  1.  The  seguence  of  events  for 
the  system  is  as  follows:  The  "slow"  field  coils 
are  energized  first,  producing  the  quasi-static 
magnetic  field  configuration  shown  in  Fig.  2. 
(Note  that  the  space  just  in  front  of  the  fast 
coil  is  a  relatively  weak  field  region  before  the 
fast  coil  is  pulsed).  Next  the  puff  valve  on  the 
diode  axis  is  suddenly  opened,  and  an  annular  gas 
puff  is  delivered  to  the  volume  in  front  of  the 
fast  coil  (at  15  cm  radius)  by  a  supersonic 
nozzle.  The  axial  puff  profile  is  sufficiently 
sharp  that  the  gas  pressure  in  front  of  the 
middle  of  the  fast  coil  can  be  up  to  100  nfTorr 
while  the  pressure  in  the  ion  diode  accelerating 
gap  4  cm  away  (to  the  right  in  Fig.  1)  is  belcw 
0.5  mTorr.  At  the  optimum  moment  relative  to  the 
time  of  opening  of  the  puff  valve  (a  function  of 
gas  species),  the  preionizer  is  energized, 
followed  about  a  microsecond  later  ty  the  fast 
coll.  The  latter  generates  a  loop  volt¬ 
age  of  17  kV  (typical)  which  drives  a  current  in 
the  preionized  gas  cloud  and  rapidly  breaks  it 
down.  The  JxB  force  on  the  plasma  current  then 
drives  the  plasna  toward  the  ion  diode  accelerat- 
ting  gap.  The  plasna  stagnates  against  the 
magnetic  field  produced  ty  the  slow  coll3,  which 
field  serves  to  magnetically  insulate  the 
ar-r-oioratlng  gap  of  the  IONGSHOT  ion  diode.  The 
precise  position  of  stagnation  is  determined  by 
the  restive  sizes  of  the  fast  and  slew  iwqnetic 
fields  and  is  designed  to  be  near  the  axial  posi¬ 
tion  of  the  metal  anode.  It  takes  about  300  ns 
for  the  plasma  to  be  driven  to  the  stagnation 


point  under  typical  plasma  source  conditions. 

The  high  voltage  pulse  is  then  delivered  to  the 
accelerating  gap. 

Figure  3  shows  the  plasna  flux  observed  on 
the  source  side  of  the  stagnation  point,  as  a 
function  of  time  between  pulsing  the  preionizer 
and  the  fast  coil.  By  selecting  this  time  in  the 
range  of  about  0.2-1.0us,  it  is  possible  to 
substantially  vary  the  arount  of  plasma  delivered 
to  the  accelerating  gap.  Varying  the  gas  puff 
pressure  in  the  range  30-70  rtfTorr,  or  the  loop 
voltage  over  the  range  14-23kV  changed  the  plasma 
flux  roughly  linearly.  The  data  in  Figure  3  were 
obtained  with  biased  Faraday  cups.  Double 
Langmuir  probe  measurements  as  a  function  of 
axial  position  across  the  stagnation  point  showed 
a  factor  of  ten  decrease  over  a  distance  of  3rtm, 
fron  a  peak  density  of  -lO'Vcm3. 

Ion  Diode  Operation  with  the  Gas-Breakdown  Plasma 
Anode 

Ion  beam  production  by  the  LONGSHOT  diode  using 
the  gas-breakdown  plasna  anode  was  monitored  with 
arrays  of  single  small-aperature  biased  (-200V) 
Farday  cups,  with  a  multiple-small-aperture  12 
cm2  biased  Faraday  cup,  with  a  4cm2  calorimeter, 
with  damage  targets  (both  with  and  without 
shadcwplates),  and  with  a  shadcwplate-Pilot-B- 
Scintillator-streak  camera  combination  to  obtain 
time  resolved  beam  optics  information,  we  first 
consider  the  basic  diode  performance  ( inpedance 
and  output  current  density)  obtained  with  the 
gas-breakdown  anode  plaana  source,  followed  by 
results  on  the  beam  quality. 

Diode  Characteristics 

By  adjustment  of  the  amount  of  plasma 
delivered  to  the  diode  accelerating  gap,  the  time 
of  arrival  of  the  plasma  relative  to  the 
application  of  the  high  voltage  pulse  from  the 
LONGSHOT  generator,  and  the  magnitude  of  the 
insulating  magnetic  field,  a  variety  of  diode 
operating  condi tons  could  be  achieved.  Figures 
4,5  and  6  shew  sets  of  data  for  pulses  which 
illustrate  the  range  of  results  that  can  be 
obtained.  Each  set  contains  inductively 
corrected  diode  voltage,  total  diode  current,  and 
two  or  more  Faraday  cup  traces.  They  show  that 
flat  or  even  rising  voltage  waveforms  can  be 
produced  by  this  ion  bean  source,  in  contrast  to 
the  typical  falling  voltage  waveform  obtained 


using  a  surface  flashover  anode  plasma  ion 
source.  Ion  current  extraction  can  be  made  to 
begin  coincident  with  the  arrival  of  the  high 
voltage  pulse  as  shown  in  Figure  5,  or  delayed 
relative  to  it  as  in  Fiqure  4b  and  4c.  Ion  beam 
pulses  of  0.6  -  lgs  above  70keV  can  be  routinely 
obtained.  In  fact,  tlie  voltage  pulses  sliown  in 
Figures  4b  and  6  are  as  long  as  "open-circuit" 
shots  without  injected  plasma.  Figure  6  shews 
that  constant  diode  impedance  could  be  maintained 
for  long  times  with  high  ion  output. 

Beam  uniformity  better  than  t  35%  in  both 
radial  and  azimuthal  directions  was  obtained  with 
Faraday  cups.  Faraday  cups  1  cm  apart  agreed  to 
♦  10%.  Gas-puff  and  preionizer  uniformity  were 
critical  to  the  attainment  of  tl  performance. 

By  averaging  Faraday  cup  traces  such  as 
those  in  Figure  5,  integrating  the  average  over 
time  and  multiplying  by  300cm?  (out  of  the  total 
400  cm7  to  correct  for  radial  profile),  we 
estimate  that  (1.5  t  0.5)x  10 17  ions  with 
energies  greater  than  60  keV  were  produced.  A  4 
cm7  calorimeter  measured  as  much  as  5J/cm2  in  the 
most  energetic  shots;  this  is  consistent  with  100 
A/cm2  for  0.5  gs,  nr  1  x  1017  lOOkeV  ions  over 
300  -  400  cm2.  Since  Faraday  cups  are  likely  to 
over-estimate  the  current  density  while  the 
calorimeter-determined  energy  density  is  a  lewer 
limit  because  of  surface  blowoff  from  the  high 
ion  fluence,  these  numbers  are  consistent.  It  is 
then  possible  to  estimate  the  ion  current 
efficiency  (ion  current/total  diode  current), 
without  correcting  for  possible  Faraday  cup 
error,  ion  current  efficiencies  up  to  70%  for  0.5 
gs  or  more  have  been  obtained.  A  lorer  limit  of 
30%  average  efficiency  . _r  the  whole  pulse  is 
given  by  the  calorimeter. 

Beam  Optics 

If  plasma  source  conditions  were  set  such 
that  the  extracted  ion  current  density,  ji  was 
<  60  A/cm2,  or  if  the  high  voltage  pulse  was 
applied  to  the  accelerating  gap  before  any  plasna 
reached  it,  streak  photographs  of  beamlets 
passing  through  an  aperture  plate  and  inpinging 
upon  a  Pllot-B  scintillator  showed  beam 
divergence  of  <3°  agreeing  with  time  integrated 
shadowbox  targets,  and  no  time  dependent  beam 
aiming  error.  However,  if  >  60  A/cm2  at  some 
time  in  some  part  of  the  diode,  and  plasma  had 
been  driven  into  the  accelerating  gap  ty  the  time 


tiie  voltaqe  pulse  arrived,  substantial  time- 
dependent  beamlet  aiming  error  occured  at  that 
location.  The  streaks  then  showed  beamlet  motion 
of  20-30°  on  a  time  scale  of  tens  to  hundreds  of 
nanoseconds.  Such  aiming  errors  might  be  due  to 
a  time-varying  uneven  anode  plasma  surface.  This 
behavior  disappeared  when  screen  mesh  was 
attached  to  the  anode  contact  forming  a  flat 
equipotential  surface  over  the  entire  anode 
area. 

Finally,  we  note  two  additional 
observations.  Because  the  puff  valve  can,  in 
principle,  be  filled  with  any  gas,  it  should  be 
possible  to  produce  pure  beams  of  many  different 
ion  species.  In  fact,  we  have  produced  more  than 
100  A/cm2  beams  using  a  nitrogen  gas  puff.  No 
major  diode  conponents  required  replacement  in 
over  700  shots  with  the  gas-breakdown  anode 
plasma  source,  confirming  the  possibility  that 
this  ion  beam  source  is  capable  of  reptitive 
pulse,  long  lived  operation. 
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Fig.  1  Diode  and  plasma  source  components:  (1) 
outer  cathode,  (2)  Inner  cathode,  (3)  Anode 
contact  and  outer  field  excluder  for  fast  field, 
(4)  Inner  field  excluder  for  fast  field,  (5) 
Outer  diode  insulating  coil,  (6)  Inner  diode 
insulating  coil,  (7)  Back  slow  coil,  (8)  fast 
coil,  (9)  Back  field  excluder  for  fast  field, 

(10)  Diode  support,  (11)  Preionizer  support,  (12) 
Mylar  collar  shield,  (13)  Stainless  steel,  (14) 
Front  nozzle  structure,  (15)  Fast  coil  insulator, 


(16)  Puff  value. 


Fig.  2. 


(b)  and  (c),  traces  for  tv*>  shots  with  the  plasma  Fig.  5  Traces  for  a  shot  with  plasma 

anode,  with  different  injection  and  timing  conditions.  injection  timing  giving  prompt  ion  beam 


t urn-on. 


Fig.  6  Diode  impedance  versus  time  for  the  shot 
with  traces  shewn. 


Fig.  3  Flux  of  plasma  from  fast  coll  as  a 
function  at  delay  between  preionizer  and  fast 
coil  firing. 
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Abstract 

The  application  of  intense  ion  rings  to  the 
heating,  stabilization  and  refluxing  of  Conpact 
Toroid  Plasmas  is  being  addressed  by  tw  experi¬ 
ments,  I REX  and  LONGSHOT.  In  IREX,  ion  rings 
containing  up  to  3  x  10 15  400keV  protons  have 
been  trapped  in  an  8  kG  magnetic  mirror.  Trap¬ 
ping  is  accomplished  by  the  use  of  a  fast  <1  us 
rise  time)  upstream  mirror  coil,  and  static  fill 
or  puffed  hydrogen  gas  to  allow  prcnpt  space 
charge  neutralization.  Trapped  rings  decay  in 
about  5  us  principally  due  to  a  combination  of 
energy  loss  and  charge  exchange  on  the  background 
gas.  In  LONGSHOT,  a  rotating  150  keV  proton  beam 
containing  up  to  10 17  particles  is  injected  into 
a  2.5  m  long  preionized  hydrogen  plasma  column  in 
a  >4.5  kG  magnetic  field  to  study  the  interaction 
of  the  ion  ring  with  the  plasma.  Results  from 
fhes<’  tvn  experiments  are  presented. 


Introduction 

The  potential  role  of  intense  ion  rings  in 
magnetic  fusion,  for  example,  to  heat,  stabilize 
and/or  reflux  Compact  Toroid  Plasna  configura¬ 
tions,  is  being  addressed  by  two  experiments  at 
Cornell,  IREX  and  LONGSHOT.  More  specifically, 
these  experiments  address  the  formation  of  ion 
rings  by  passing  an  ion  beam  from  an  annular 
magnetically  insulated  intense  ion  beam  diode 
through  a  magnetic  cusp,  the  trapping  of  such 
rings  in  a  magnetic  mirror  configuration  (IREX), 
and  the  merging  of  these  rings  with  a  Compact 
Toroid  Plasma  (LONGSHOT) .  This  paper  discusses 
the  progress  made  in  these  two  experiments  since 
the  last  conference  in  this  series  (1).  In  IREX, 
we  have  trapped  up  to  3  x  10 15  400keV  protons  in 
an  8  kG  magnetic  mirror.  In  LONGSHOT,  we  are 
studying  the  propagation  of  150keV  ion  rings  in  a 
>.  1/2  m  long  pre ionized  hydrogen  plasma  in 
preparation  for  an  experiment  to  merge  a  proton 
ring  with  a  compact  Tbroid  Plasma. 


IREX 

The  Ion  Ring  Experiment  (IREX)  is  shown 
schematically  in  its  present  form  in  Fig.  1.  The 
annular  magnetically  insulated  diode  shown  in 
Fig.  I  is  used  to  inject  a  >430  keV  proton  beam 
thrcugh  a  magnetic  cusp  produced  by  the  super¬ 
position  of  the  field  f ran  the  solenoid  sur¬ 
rounding  the  experiment,  and  the  diode  insulating 
field  produced  by  the  coil  in  the  anode.  The 
resulting  rotating  proton  beam  ("ring")  has  80- 
90%  of  its  energy  in  the  transverse  direction. 
Trapping  i3  accomplished  by  injecting  the  axis- 
encircling  proton  ring  through  a  fast  risetime 
!lys)  upstream  mirror  coil  capable  of  achieving  a 
1.4:1  peak  mirror  ratio,  and  reflecting  the  ring 
off  a  static  1.23-1.45:1  mirror  located  1.5m 
downstream.  In  order  to  provide  space  charge 
neutralization  for  the  ions,  the  rings  are  pro¬ 
pagated  and  trapped  in  a  static  >50nm>rr  H2  gas 
fill,  or  in  gas  puffed  only  in  regions  of  spati¬ 
ally  varying  magnetic  field.  The  <10-U  Torr  ion 
diode  region  is  separated  from  the  gas  filled 
region  of  the  experiment  by  a  2i*n  mylar  foil. 

The  protons  are  detected  using  magnetic  probes, 
Faraday  cups,  heat  sensitive  "witness  targets", 
and  CR-39  particle  track  detector  material . 

Initial  Trapping  Experiments 

The  earliest  trapping  experiments  (2)  were 
performed  with  1.23:1  static  downstream  mirror,  a 
1.22:1  peak  pulsed  upstream  mirror,  and  a 
relatively  weak  but  reproducible  proton  ring 
containing  typically  2xl015  protons  on  its  first 
pass  through  the  mirror.  The  pulsed  power  gen¬ 
erator  driving  the  ion  diode  was  the  7n,  100  ns. 
Neptune  generator  we  have  previously  described 
|3J.  The  observed  diamagnetism  of  the  injected 
rings  in  the  middle  of  the  mirror  (z  *  136  cm) 

was  up  to  4B/B  of  about  0.01,  where  B  ,  the 
o  o 

applied  magnetic  field,  was  typically  8  kG. 

In  general,  hydrogen  fill  pressures  above 
100  rrfrbrr  gave  good  propagation,  reflection  and 
trapping,  while  fill  pressures  below  50  nirorr 


gave  poor  results.  Good  trapping  was  character¬ 
ized  'v  detectable  ion  rings  (>0.5A/cm2)  for  6-8 
passes  of  the  ring  through  the  length  of  the 
mirror,  a  time  of  about  4  us  (about  50  times  the 
ion  cyclotron  period).  The  estimated  initial 
number  of  trapped  protons  (i.e.  -  reflected  off 
the  downstream  and  then  the  pulsed  upstream 
mirrors)  was  typically  4  x  lO14,  decreasing  to 
2  x  10ls  on  the  eighth  pass.  The  particle  loss 
is  believed  to  have  been  due  to  a  combination  of 
charge  exchange  and  mirror  losses,  the  latter 
because  of  the  oscillatory  upstream  mirror  (see 
Fig.  3,  upper  trace).  A  data  corrposite  for  a 
typical  ion  ring  trapped  in  a  static  hydrogen 
fill  was  shown  in  Ref.  2.  Here  (Fig.  2)  we 
present  a  data  composite  for  a  trapped  ring  in 
which  localized  gas  clouds  were  produced  in  the 
JRFX  experiment  chamber  by  three  puff  valves, 
tv->  located  in  the  transverse  magnetic  field 
region  near  the  injector  and  one  just  outside  the 
downstream  magnetic  mirror.  The  H2  gas  pressure 
profile  for  the  trapped  ring  shown  in  Fig.  3  is 
presented  in  Fig  4. 

Puffed  gas  experiments  were  performed  with 
two  purposes  in  mind:  1.)  to  determine  what 
reqions  of  the  experiment  required  gas,  and  at 
what  pressure,  to  trap  a  ring,  and  2).  to  in¬ 
crease  the  ring  lifetime  with  respect  to  charge 
exchange  losses  by  eliminating  the  gas  where  it 
was  not  needed.  The  Faraday  cup  traces  in  Fig.  3 
shew  a  detectable  ion  ring  up  to  nearly  4uS, 
about  the  same  duration  as  with  a  >  100  irfTorr 
static  fill  during  these  initial  trapping  experi¬ 
ments.  When  the  peak  of  the  puffed-gas  profile 
in  the  cusp  region  was  >75  nfTorr,  there  was  only 
modest  degradation  of  the  ring  injected  into  the 
TREX  mirror  region,  as  determined  by  a  decrease 
in  the  muter  of  injected  protons  and  by  addi¬ 
tional  spreading  of  the  ring.  When  the  gas  puff 
peak  was  reduced  to  40-75  irtTorr,  the  injected 
ring  peak  diamagnetism  at  z  *  136  an  was  reduced 
by  a  factor  of  three,  its  width  increased  by  50*, 
and  the  estimated  number  of  particles  decreased 
by  a  factor  of  four,  relative  to  the  ring  shown 
in  Fig.  2.  It  was  determined  that  gas  was  needed 
throughout  the  region  of  nonuniform  magnetic 
field  at  the  Injector  end  of  IREX.  The  strength 
of  the  rings  reflected  fran  the  downstream  mirror 
(as  determined  ty  the  peak  current  density  of  the 
reflected  rlnq  at  z  =  136  cm)  was  also  affected 
by  the  gas  puff  pressure,  but  only  up  to 
15-20  mtorr,  above  which  no  further  increase  in 


ring  strength  was  obtained.  Below  15  mTot r,  the 
reflected  ring  current  density  was  roughly 
linearly  dependent  on  puff  gas  pressure.  No 
increase  in  trapped  ring  lifetime  was  seen  in 
these  experiments.  This  is  probably  because  of 
the  poorly  crowbarred  upstream  mirror  in  these 
initial  experiments  (see  upper  trace  in  fig.  2), 
and  the  fact  that  scattering  and  charge  exchange 
probably  were  not  decreased  by  more  than  a  factor 
of  2  -  3  using  the  puff  gas  profile  shown  in 
Fiq.  3. 

Experime nts_  with  Stronger  Injected  Rings 

The  injected  rings  used  for  the  initial 
trapping  experiments  were  optimized  for  repro¬ 
ducibility.  An  additional  series  of  experiments 
[41  was  performed  in  which  the  number  of  in¬ 
jected  protons  was  optimized  within  the  capa¬ 
bilities  of  the  7  n  pulsed  power  generator  then 
being  used.  In  addition,  the  downstream  magnetic 
mirror  ratio  was  boosted  to  a  neximum  of  1.45:1, 
and  the  pulsed  upstream  mirror  was  improved  by 
both  a  greater  peak  mirror  ratio,  1.4:1,  and  a 
better  crowbar.  This  revised  magnetic  field 
configuration,  which  is  the  one  shewn  in  Fig.  1, 
has  been  used  only  with  static  hydrogen  fills. 
Injected  rings  of  up  to  SB/Bq  =  0.04  were 
obtained  at  z  =>  136  an.  With  the  1:45.1  down¬ 
stream  mirror,  up  to  90%  of  the  protons  were 
reflected,  producing  a  reflected  ring  4B/Bq  of 
more  than  0.02  at  the  same  central  location. 

Since  the  magnetic  probes  that  measure  6B  are 
shielded  from  the  ring  diamagnetism  after  two  or 
three  bounces  by  currents  induced  in  the  beam- 
formed  plasma,  we  have  had  to  rely  on  magnetic¬ 
ally  insulated,  biased  (-200V)  Faraiday  cups  to 
monitor  the  strength  of  trapped  rings.  Figure  4 
shows  sanple  data  obtained  from  a  trapped  ring. 
Note  that  the  Faraday  cup  is  collimated  and  aimed 
to  view  only  dcwnstrean  moving  rings.  At  least 
seven  distinct  ring  current  density  peaks  are 
visible  on  the  signal,  indicating  that  many  round 
trip  passes  through  the  mirror,  or  fourteen 
mirror  reflections,  during  its  lifetime.  Notice 
that  the  magnetic  probe  signal  peaks  drop  off 
more  quickly  than  the  Faraday  cup  peaks.  The 
peak  current  density  of  the  entering  ring,  cut 
off  in  the  trace  shown  was  about  20A/cm2 .  Notice 
that  the  magnetic  probe  singal  peaks  drop  off 
more  quickly  than  the  Faraday  cup  peaks.  In  the 
better  shots  in  this  series,  3x10 15  or  more 
protons  (J200.7)  were  estimated  to  be  trapped  in 


the  magnetic  mirror,  with  life  times  (using  75  - 
100  rrtTorr  of  H2)  of  4-5  us.  Two  more  conplete 
data  sets  of  a  trapped  ring  frcm  which  such 
numbers  can  be  obtained  are  shown  in  ref. 4.  The 
ensemble  of  such  trapped  ring  data  supports  the 
hypothesis  that  once  the  ring  is  trapped  (i.e.- 
after  the  first  bounce  off  the  upstream  mirror), 
charge  exchange  losses  are  responsible  for  the 
decay  of  the  ring.  (See  ref.  4  for  details). 
Finally  notice  that  Fig.  4  suggests  that  in  addi¬ 
tion  to  the  clear  successive  peaks  of  the  ring 
bouncing  in  the  mirror,  there  is  a  uniform  few 
A/cm5  "background"  ion  ring  filling  the  entire 
mirror  after  the  second  bounce.  These  two  com¬ 
ponents  appear  to  decay  at  the  same  rate,  as  they 
should  if  charge  exchange  is  the  loss  mechanism. 

Fxperiments  with  an  Upgraded  Generator _ 

A  third  trapping  run  with  the  ion  diode 
driven  by  a  30  version  of  the  Neptune  generator 
is  presently  in  progress  on  IREX  using  static 
hydrogen  gas  fills.  Although  the  upstream  mirror 
has  not  yet  been  used,  we  have  seen  current 
densities  of  injected  rings  and  rings  reflected 
off  the  1 .45 s 1  downstream  mirror  at  least  twice 
as  large  as  in  the  previous  run,  as  illustrated 
by  the  Faraday  cup  traces  shewn  in  Fig.  5.  The 
peak  current  density  shown  inplies  a  ring  of  peak 
*B/B0  of  order  0.1,  and  is  consistent  with  dia¬ 
magnetic  signals.  This  injected  ring  is 
reflected  from  the  downstream  mirror  as  well  as 
weaker  rings  vrere,  so  we  expect  soon  to  be  able 
to  study  trapped  rings  of  at  least  double  the 
previous  current  and  energy. 

LONGS HOT 

The  LONGSHOT  II  proton  ring  -  plasma 
interaction  experiment  is  shown  schematically  in 
Fig.  6.  The  LONGSHOT  II  Injector  with  a 
preionized  anode  plasma  ion  source  (described  in 
detail  in  a  paper  fcy  Greenly,  Ueda,  et  ai.  in 
this  conference)  produces  an  annular  beam 
typically  containing  more  than  .10* 7  >60keV 
protons  (>  1  kJ).  This  beam  is  injected  through 
a  cusp-iike  magnetic  field  into  a  2.5m  long,  15 
cm  diameter,  preionized  hydrogen  plasma,  produced 
fcy  a  Z-dlscharge  in  the  presence  of  >4.5kG 
solenotdal  magnetic  field.  To  form  the  plasma,  a 
damped  oscillatory  circuit  (C*100pF,  L”3.2  yH,  R 
”  Sftnfl)  drives  current  between  annular  stainless 
stoel  electrodes  separated  ty  about  2m  In  a 


puffed  neutral  hydrogen  fill  of  about  5  sfTorr  for 
one  half  cycle  (65  us),  after  which  it  is  crow- 
barred.  The  afterglow  plasma  Into  which  the  ion 
ring  is  injected  has  been  determined  by  double 
Langmuir  probes  and  a  microwave  interferometer  to 
be  a  peak  density  of  10  ’Von3  and  a  temperature 
of  about  3eV.  An  80%  transparent,  fine  stainless 
steel  mesh  stretched  across  the  entrance  to  the 
interaction  region  (at  z=20cm)  and  "grounded"  to 
the  drift  chamber  walls  prevents  injector 
operation  from  being  affected  fcy  plasma  drifting 
into  it  from  the  Z-dlscharge. 

Experiments  with  this  ring-plasma  have 
recently  begun.  Early  results  indicate  that  the 
ring  is  entering  the  plasma,  and  the  diamagnetism 
of  the  ring  can  be  seen  on  axis  in  the  plasma 
near  the  injector.  The  magnetic  signals  decrease 
rapidly  In  lm  of  axial  distance  through  the 
plasma.  The  magnetic  signal  outside  the  ring 
radius  at  the  vacuum  wall  shows  strong  evidence 
of  induced  plasma  current,  and  oscillatory 
signals  persist  after  passaqe  of  the  ring. 

These  experiments  will  study  the  ring-plasma 
interaction,  looking  in  particular  for  evidence 
of  slowing  of  the  ring  axial  velocity  by  coupling 
ring  energy  to  magnetosonic  waves  in  the 

plasma. (5) 
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capacitor  supplying  (C)  pulsed  upstrean  mirror; 
(D)  power  supply  for  (E)  solenidal  field;  ( F ) 
supplemental  downstrean  mirror  coil;  (G)  typical 
Faraday  cup  location; (L)Neptune  generator 
supplying  ion  diode  anode  and  (N)  cathode  with 
2um  foil.  Net  applied  magnetic  field  on  axis  is 
plotted  below.  Diode  detail  is  also  shown. 
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Fig.  2  Data  for  an  IREX  shot  into  puffed  gas. 

All  traces  are  500ns/div;  the  upstream  gate 
mirror  strength  is  shown.  Other  traces  are  ion 
ring  diamagnetic  signal  N  (shaded,  below  the 
abcissas)  at  axial  locations  indicated  at  right. 
The  two  traces  above  the  absissa  the  the  136  cm 
position  are  signals  from  Faraday  cups  collecting 
(solid  line)  all  ring  ions,  and  (dashed  line) 
only  upstream-moving  (reflected)  ring  ions. 
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Fig.  4  ion  ring  trapping  with  7  fl  injector  pulse 
shown  by  diode  voltage  and  current  traces.  Ion 
current  density  at  160  an  axial  position, 
measured  by  Faraday  cup,  and  ring  diamagnetic 
signal  at  the  vacuum  wall,  showing  rapid 
shielding  of  the  ring  self-field  by  induced 
currents  in  the  ring-formed  background  plasma. 
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Fig.  3  Puffed-gas  pressure  versus 
axial  position  for  shot  of  Figure  2 
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Fig.  6  LONGSHOT  II-  Proton  Ring 
Plasma  Interaction  Experiment : 

(A)  to  generator;  (B)  ion  diode; 

(C)  grounded  grid;  (D)  half-cusp 
producing  alumimm  plate;  (E)  protor 
orbit  coup  re ss ion  coils;  (F) 
solenodial  field  producing  coils; 

(G)  ground  Z -discharge  electrode; 

(H)  high-voltage  z-dlscharge 
electrode;  (I)  puff  valve  for  Z-dis- 
charge  fill. 
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propagation  OF  INTENSE  ION  BEAMS  IN 
STRONG  MAGNETIC  FIELDS  AND  BACKGROUND  PLASMA 

R.  Kraft*,  B.R.  Kusse  and  J.J.  Moschella 
Cornell  University,  Ithaca,  New  York 


Abstract 

The  dynamics  and  propagation  of  an  intense 
ion  beam  under  varing  magnetic  field  and  back¬ 
ground  plasma  conditions  have  been  studied  in  a 
series  of  experiments.  The  70  nanosecond  (fwhm) 
beam  consisted  of  8  KA  of  360  KeV  protons.  The 
first  set  of  experiments  investigated  theoretical 
models  for  ion  beam  propagation  and  neutralization. 
Conventional  diagnostics  and  a  new  method  for  in¬ 
terpreting  biased  charge  collector  signals  were 
used.  The  second  study  concerned  the  focusing  of 
an  intense  ion  beam  as  it  propagated  from  a  field 
free  region  into  a  strong  magnetic  solenoid.  The 
third  set  of  experimants  concentrated  on  the  ob¬ 
servation  of  microwaves  produced  by  ion  beam- 
plasma  interactions. 

Experimantal  Apparatus 

The  experimantal  apparatus  is  shown  in  Fig.  1. 
The  Ion  beam  was  produced  by  a  planar,  race  track 
type,  magnetically  insulated  diode  driven  by  a 
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marx  generator  and  pulse  forming  network.  The  pro¬ 
ton  beam  had  a  10  cm  x  13  cm  cross-section  which 
could  be  collimated  down  to  a  4  cm  diameter  by  two 
aperture  plates.  After  collimation  the  beam  enter¬ 
ed  a  1  m  drift  region  where  axial  and  transverse 
magnetic  fields  could  be  applied.  A  conical  9- 
pinch  plasma  gun  was  located  at  the  downstream  end 
of  this  drift  region  and  was  used  to  provide  a  pre¬ 
ionized  medium  into  which  the  beam  could  be  propa¬ 
gated.  A  puff  gas  valve  located  in  the  plasma  gun 
could  be  fired  without  voltage  on  the  0-plnch  coil 
to  establish  an  adjustable,  background  neutral 
pressure  In  the  drift  region. 

Beam  Propagation  and  Neutralization 
Theoretical  Picture  It  has  been  proposed 


that  the  ion  beam,  as  it  Initially  leaves  the 
cathode,  picks  up  a  neutralizing  electron  distri¬ 
bution  which  is  composed  of  two  delta  functions, 
one  at  zero  velocity  and  one  at  twice  the  ion 
velocity  [1].  The  total  number  of  electrons 
equals  the  total  number  of  ions  and  the  beam  is 
current  neutralized.  It  is  also  essentially 
space  charge  neutralized  except  for  the  over¬ 
shooting  electrons  which  establish  a  self-consis¬ 
tent  positive  beam  space  charge  potential. 

This  electron  distribution  is  unstable  and 
is  predicted  to  quickly  thermal ize  into  a  drift¬ 
ing  Maxwelllam  centered  at  the  ion  velocity  [2,3] 
This  thermal izat ion  also  lowers  the  beam  space- 
charge  potential. 

Biased  Charge  Collector  When  used  with  in¬ 
tense  ion  beams,  the  aperature  of  a  biased  charge 
collector  reduces  the  ion  current  drematically 
but  not  enough  to  allow  one  to  Ignore  space 
charge  effects  In  the  region  between  the  apera¬ 
ture  and  collector.  Including  these  effects  in 
an  analysis,  it  is  possible  to  develope  a  set  of 
I-V  curves  for  a  collector  which  allows  determin¬ 
ation  of  the  neutralizing  electron  temperature  as 
well  as  the  ion  beam  density. 

The  geometry  of  the  collector  is  shown  in 
Fig.  2  and  consists  of  aperature  and  collector 
plates  and  a  floating  bias  supply.  The  current 


signal  is  measured  as  a  voltage  across  the  50  ohm 
resistor.  The  Incident  beam  Is  assumed  to  be  com 
posed  monoenergetic  ions  at  density  NQ  and  vel¬ 
ocity  vQ  together  with  the  neutralizing  electrons 
of  equal  density  but  characterized  by  a  drifting 
Maxwellian  centered  at  vQ  and  a  temperature  Tg. 
The  analysis  for  the  I-V  characteristic  for  the 
collector  proceeds  from  assuming  that  the  apera¬ 
ture  plate  attains  the  floating  potential  with 


respect  to  the  beam  potential  such  that  the  in¬ 
cident  ion  and  electron  currents  are  equal.  A 
fixed  voltage,  VQ,  is  applied  between  the  aperature 
and  the  collector  plates.  The  ions  are  assumed  to 
be  unperturbed  and  their  current  density  and  vel¬ 
ocity  are  taken  to  be  constant  both  inside  and 
outside  the  collector.  A  space  charge  potential, 
developes  in  the  column  between  the  aperature 
and  collector  plates  to  keep  the  electron  and  ion 
densities  equal.  Once  *)fc  is  known  the  electron 
distribution  and  current  density  at  the  collector 
plate  can  be  determined  [4], 

In  obtaining  the  electron  distribution  two 
conditions  must  be  considered.  If  the  collector 
potential  dc=  df-  VQ  is  greater  than  the  ion  col¬ 
umn  potential  there  are  no  reflected  electrons 
at  the  collector  and, 
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must  be  added  to  the  electron  density.  Here  fg(v) 
is  the  same  function  given  by  Eq.  1.  For  the  case 
with  *i1c>  j)c  the  column  potential  is  again  deter¬ 
mined  by  requiring  the  electron  density  to  equal 
Nq.  Once  is  known,  the  electron  current  and 
therefore  the  net  current  to  the  collector  can  be 
calculated  as  a  function  of  VQ,  vQ,  Te  and  NQ. 

A  set  of  curves  showing  the  collector  current 
as  a  function  of  the  applied  voltage  is  shown  in 
Fig.  3  where  Tg  has  been  used  as  the  parameter  to 


Mg.  3  Noma  1< z«d  Collector  Current  »i.  tin  Voltage 

generate  the  set.  It  can  be  seen  that  the  curves 
behave  as  expected  until  V0  reaches  a  critical 
value  below  which  the  space  charge  effects  take 


over  and  clamp  the  current.  This  critical  potent¬ 
ial  depends  on  Tg  and  consequently  curve  fitting 
experimantal  data  allows  a  determination  of  Tg 

as  well  as  N  . 

o 

Measurements  of  and  Tg  A  capacitive 
probe  shown  in  Fig.  4  and  a  biased  charge  col¬ 
lector  were  used  to  measure  T  ,  N  and  d.  .  The 

e  o  b 


riq.  4  Capacitive  Probe 


capacitive  probe  was  designed  to  sit  at  the  float¬ 
ing  potential  and  was  obtained  from  d^  using 

T  . 
e 

The  experimental  configuration  was  as  shown 
in  Fig.  1  with  a  permanent  magnet  and  pole  pieces 
to  establish  a  localized,  transverse  field  of  ap¬ 
proximately  200  G.  Immediately  down  stream  from 
this  transverse  field  was  a  grounded  conducting 
screen  which  could  supply  a  new  set  of  neutraliz¬ 
ing  electrons  to  replace  those  turned  by  the 

field.  Potential  and  T  measurements  were  made 
e 

down  stream  from  this  screen  and  a  plot  of  6^  vs. 
distance  from  the  screen  Is  shown  In  Fig.  5.  The 
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electron  temperature,  after  thermal izat ion,  was 
measured  to  be  156  eV  and  the  predicted  decrease 
of  the  potential  with  distance  is  clear. 

Collective  Ion  Beam  Focusing 

Because  of  the  intensity  of  the  neutralized 
ion  beam  considered  in  this  work,  magnetic  focus¬ 
ing  can  require  collective  effects  to  adequately 
describe  the  phenomenon.  Collective  focusing 
using  a  short  magnetic  lens,  where  the  focal  point 
occurred  in  a  field  free  region  down  stream  from 
the  lens,  has  been  described  by  Robertson  [5], 


Our  studies  involved  collective  focusing  of  the 
ion  beam  which  occurred  as  the  beam  entered  a  long 
magnetic  solenoid. 

Experimental  Set-up  The  experiment  was  ar¬ 
ranged  as  shown  in  Fig.  1  with  the  solenoidal  coil 
energized.  The  neutralized  ion  beam  was  collimat¬ 
ed  and  propagated  a  short  distance  through  a  field 
free  region  before  it  encountered  the  radial 
fields  of  the  solenoid.  Within  the  solenoid  the 
biased  charge  collector  was  used  to  measure  the 
radial  and  axial  profiled  of  the  beam  intensity. 

Focusing  Theory  The  collective  focusing  ef¬ 
fect  involves  a  co-operatinn  between  the  beam  ions 
and  their  neutralizing  electrons.  Single  particle 
theories  do  not  apply  because  of  the  large  space 
charge  fields  which  can  build  up.  The  parameters 
of  the  experiment  were  such  that  the  electrons,  by 
themselves,  could  not  have  penetrated  the  radial 
fields.  The  ions  had  no  trouble  crossing  and 
their  space  charge  was  able  to  add  energy  to  the 
electrons  and  pull  them  across.  The  ions,  by 
themselves,  would  have  focused  at  a  position  In 
the  solenoid  equal  to  the  distance  they  could 
travel  in  1/2  a  cyclotron  period.  As  single  par¬ 
ticles  the  electrons  wanted  to  focus  in  a  much 
shorter  distance.  As  these  events  try  to  occur 
radial  electric  fields  built  up  resulting  in  a 
convnon,  collective  focus  inbetween  the  single  par¬ 
ticle  electron  and  ion  focal  points. 

This  process  can  be  described  analytically  by 
assuming  quasi-neutral ity  and  followinq  a  calcul¬ 
ation  similar  to  Robertson's.  The  Lagrangians  for 
the  electrons  and  ions  are  by, 

L  *  1/2  m  v2  +  q  (7  *  S)  -  q  0  (3) 

ft  a  a  a  a  a  '  ’ 


where  the  a  subscript  is  used  to  indicate  ions  or 
electrons,  m  is  the  mass,  qg=  -e,  v  the  vector 
velocity,  A  the  vector  potential  and  0  the  scalar 
potential.  Quasi -neutrality  requires  that  the  rad¬ 
ial  and  axial  motion  of  the  electrons  and  ions  be 
equal  and  the  following  equations  of  motion  result. 
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In  Fig.  6a  radial  trajectories  of  monoenerget- 
1c  ions  are  plotted.  As  long  as  the  ions  remain 
In  the  solenoid,  multiple  focusing  is  predicted. 

In  Fig.  6b  the  radial  trajectories  for  ions  with  a 
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velocity  spread  have  been  plotted.  While  the 
first  focal  point  is  still  present,  subsequent 
ones  are  smeared  out  and  indistinguishable. 

Focusing  Measurements  A  biased  charge  col¬ 
lector  was  used  to  measure  the  peak  current  den¬ 
sity  as  a  function  of  axial  position  within  the 
solenoid.  The  results  are  shown  in  Fig.  7a.  On 
the  same  graph  is  a  similar  measurement  with  no 
solenoidal  field  to  indicate  the  normal  beam  div¬ 
ergence.  Only  a  single  focal  point  was  observed. 
The  radial  profile  at  the  focal  point  is  shown  in 
Fig.  7b.  The  measured  and  theoretical  focal 
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points  as  a  function  of  magnetic  field  strength 
are  shown  in  Fig.  8. 

The  focusing  effect  was  very  dependent  on  the 
ions  being  able  to  pull  original  neutralizing  el¬ 
ectrons  all  the  was  across  the  field.  If  it  was 
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easier  for  the  ions  to  pick  up  new  electrons  as  a 
result  of  an  electron  source  in  the  radial  field 
region  (i.e.  a  screen  or  wall  bombardment  by  an 


uncollimated  ion  beam)  or  a  plasma  prefill,  the 
focus  was  destroyed. 

Ion  Beam-Plasma  Microwave  Interactions 

It  is  believed  that  microwave  Interactions 
play  a  large  role  in  ion  beam  neutralization  and 
propagation  [2,3].  In  order  to  propagate  an  in¬ 
tense  ion  beam  in  vacuum  a  large  degree  of  magne¬ 
tic  and  electrostatic  neutralization  is  required. 
This  neutralization  of  the  cold  drifting  ion  beam 
is  provided  by  an  electron  population  whose  vel¬ 
ocity  distribution  evolves  with  time.  Initially, 
as  the  electrons  are  drawn  into  the  beam  channel, 
they  are  cold  and  attain  a  drift  velocity  greater 
than  that  of  the  ions.  As  they  overshoot  the  ions 
they  are  decelerated  to  approximately  zero  veloc¬ 
ity.  The  two  component  distribution  that  results 
is  unstable  and  quickly  thermal izes  by  means  of  a 
microwave,  two-stream  instability  to  a  hot,  drift¬ 
ing  Maxwellian  centered  at  the  ion  velocity.  If 
this  neutralized  beam  then  propagates  through  a 
background  plasma  further  two-stream  interactions 
can  be  present. 

Observations  of  Instabilities  Glass  enclosed 
electrostatic  probes  were  used  to  observe  micro- 
wave  signals  both  inside  and  outside  the  beam  chan¬ 
nel  of  the  experiment  shown  in  Fig.  1.  The  strong, 
broadband  microwave  signals  from  the  diode  made 
measurements  during  the  Initial  neutralization  and 
thermal izat ion  near  the  cathode  impossible.  Our 
cleanest  measurements  occurred  when  we  stripped 
the  beam  of  its  neutralizing  electrons  and  forced 
it  to  pick  up  new  cold  ones  well  down  stream  from 
the  cathode.  Under  these  conditions  we  were  able 
to  observe  an  interaction  between  the  cold  elect¬ 
rons  and  those  in  a  background  plasma. 

This  interaction  was  strongest  when  the  beam 
was  used  to  form  a  plasma  from  a  neutral  prefill 
and  cold  fast  electrons  were  produced  as  the  beam 
entered  the  magnetic  solenoid.  With  this  back¬ 
ground  plasma  no  collective  focusing  occurred.  A 
plot  of  the  frequency  of  the  microwave  signal  as  a 
function  of  the  prefill  pressure  is  shown  In  Fig. 

9.  The  pressure  was  measured  in  the  center  of  the 
center  of  the  solenoid  at  the  time  of  the  beam 
pulse  and  can  be  converted  to  a  plasma  density 
through  a  rate  of  ionization  equation.  The  fre¬ 
quency  was  determined  by  an  x-band  dispersive  line 
and  scaled  with  the  background  pressure  and  plasma 
density  but  was  independent  of  magnetic  field 
strength.  The  microwave  signal  strength  as  a  fun¬ 
ction  of  position  is  shown  in  Fig.  10.  The  z  =  0 


Fig.  9  Microwave  Frequency  vs.  Neutral  Pressure 


position  is  under  the  first  solenoid  coil  and  the 
spatial  growth  indicated  a  convective  instability 


with  a  20  cm  e-folding  length. 

Conclusions 

Several  experiments  confirming  the  theoret¬ 
ical  picture  for  neutralized  ion  beam  propagation 
have  been  performed.  Evidence  of  the  production 
of  fast,  cold,  neutralizing  electrons  that  ther¬ 
mal  ize  through  a  two-stream  instability  has  been 
presented.  Observations  of  the  beam  potential 
show  its  expected  decrease  as  the  electrons 
thermalize  into  a  hot,  drifting  Maxwellian  dis¬ 
tribution.  This  distribution  was  measured  using 
a  new  technique  for  analysing  biased  charge  col¬ 
lector  signals.  No  microwaves  driven  by  the  neg¬ 
ative  slope  of  the  ion  distribution  were  observed. 
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Abstract 


We  have  measured  the  electric  field  distribution 
across  the  gap  of  magnetically  insulated  intense  ion  diodes 
by  observing  the  Stark  shift  of  ion  line-emission  from  the 
gap.  The  measured  distribution  yielded  the  time  depen¬ 
dence  of  the  actual  diode  gap,  which  decreased  rapidly 
early  in  the  pulse.  The  electron  density  distribution,  in¬ 
ferred  from  the  measured  E(z),  was  compared  to  one¬ 
dimensional  calculations  using  the  measured  gap.  Elec¬ 
trons  were  not  confined  to  a  well  defined  sheath  near  the 
cathode,  but  were  spread  significantly  towards  the  anode, 
especially  in  vaned  cathode  diodes.  These  results  serve  to 
explain  the  enhancement  of  the  measured  ion  current  den¬ 
sity  over  the  space-charge-limited-value  (calculated  using 
the  measured  gap)  that  we  observed.  We  also  measured 
the  velocity  distribution  parallel  to  the  magnetic  field  lines 
of  C++  and  A1++  ions  accelerated  in  the  gap  by  observ¬ 
ing  the  line  Doppler-broadening.  This  velocity  component 
was  considerably  smaller  than  the  transverse  proton  ve¬ 
locity  previously  measured  outside  the  diode. 

Introduction 

A  true  understanding  of  the  operation  of  high  power 
diodes  requires  measurements  of  their  detailed  properties, 
such  as  the  actual  accelerating  gap  spacing,  the  charge 
density  distribution,  the  properties  of  the  anode  plasma, 
etc.  However,  the  small  gap  spacings  (typically  <  lem), 
and  large  applied  electric  fields  (typically  >  lA/V/em) 
make  direct  measurements  within  high  voltage  gaps  very 
difficult.  In  this  paper  we  report  spectroscopic  measure¬ 
ments  of  the  electric  field  as  a  function  of  distance  from 
the  anode,  and  of  the  ion  transverse  velocities  within  a 
magnetically  insulated  ion  diode.  These  measurements 
were  made  by  observing  spontaneous  line-emission  from 
ions  excited  in  the  anode  plasma  and  accelerated  in  the 
diode  gap.  The  electric  field  measurements  reveal  the 
actual  accelerating  gap  spacing  by  locating  the  bound¬ 
aries  of  the  electric-field-excluding  anode  and  cathode 
plasmas,  rather  than  the  position  of  the  ~  10,8em~3  or 
higher  plasma  boundaries  observed  in  previous  studies  (1- 
3).  These  measurements  also  give  information  about  the 
electron  density  within  the  diode,  which  was  found  to  be 
inconsistent  with  theoretical  predictions  assuming  a  well 
confined  electron  sheath  near  the  cathode  (4).  The  mea- 
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sured  gap  spacing  and  electron  density  distribution  serve 
to  explain  the  enhancements  we  have  measured  in  the  ion 
current  densities  over  calculated  [4]  values. 

The  measurements  reported  here  were  initiated  as  a 
result  of  a  recent  suggestion  [5|  that  the  electric  field  in 
the  acceleration  gap  of  a  high  power  diode  can  be  deter¬ 
mined  by  measuring  the  Stark  shift  of  line-emission  from 
especially  selected  ions  in  the  gap,  as  illustrated  in  Fig. 
la.  Observations  at  different  distances  z  from  the  an¬ 
ode  yields  the  electric  field  distribution  E(z)  across  the 
gap.  With  this  method,  the  measurement  integrates  over 
the  field  of  view  in  front  of  the  anode.  The  line-emission 
(collected  parallel  to  the  electrodes)  is  also  Doppler  broad¬ 
ened  due  to  the  transverse  ion  motion  (parallel  to  the  elec¬ 
trodes).  Therefore,  as  discussed  in  (5),  the  measurement 

of  the  ion  transverse-velocity  distribution  in  the  diode  gap 
is  accomplished  by  observing  the  spectral  profile  of  a  line 
which  is  negligibly  shifted  by  the  electric  field  relative  to 
the  Doppler  broadening.  We  selected  the  Al++4rfs/j  -♦ 
ips/2  emission  (4529  A)  for  the  electric-field  measure¬ 
ments,  and  confirmed  it  with  the  Al++4d3/j  — *  4pi/2  line 
(4512  A).  The  Al++4p  -  3d  (3602  A)  line  as  well  as  the 
C++  2s(*S)3p  — *  2s3 a  (4647  A)  line  were  used  for  the 
transverse  velocity  measurements. 

Apparatus  and  Measurements 
Experimental  arrangement 

We  used  the  planar  magnetically  insulated  diode  [6j 
which  is  shown  in  Fig.  lb.  The  dielectric-anode  surface 
was  made  up  of  a  140x50x1.6  mm*  polyethelene  sheet 
with  the  long  dimension  parallel  to  the  applied  magnetic 
field  Bt.  Line-radiation  emitted  in  the  s  direction  from  a 
planar  surface  parallel  to  the  anode  was  directed  by  the 
mirror  M  to  a  lens  L  which  focused  it  on  the  spectrometer 
input  slit.  Using  the  cylindrical  lens  CL,  the  spectrom¬ 
eter  output  was  magnified,  then  focused  on  an  array  of 
seven  fiber  bundles  placed  parallel  to  the  spectrometer 
slit.  Each  bundle  transmitted  the  light  to  a  PM  tube,  the 
trace  of  which  was  recorded  with  a  temporal  resolution 
of  8ns.  The  spatial  resolution  in  the  x-directlon  and  the 
spectral  resolution  were  0.6  mm  and  0.7  A,  respectively. 
Thus,  seven  points  of  the  spectral  line  profile  were  ob¬ 
tained,  giving  the  electric  field  or  the  transverse  velocity 
distribution  as  a  function  of  time  in  a  single  discharge. 


Fig.  1.  (a)  Method  for  measuring  the  electric  field  E 
in  the  acceleration  gap  of  a  magnetically  insulated  diode, 
(b)  Schematic  illustration  of  the  planar  magnetically  insu¬ 
lated  diode  and  the  optical  arrangement.  The  distance  of 
the  observation  region  from  the  anode  is  varied  by  moving 
the  mirror  M  in  the  x  direction. 

Measurements  of  the  electric  field 

The  profile  of  the  unshifted  4529  A  Al++  emission 
from  the  anode  plasma,  and  the  red  shifted  emission  from 
the  gap  are  shown  in  Fig.  2.  The  signal  from  the  gap 
at  the  zero-field  wavelength  (i.e.,  the  signal  in  the  fiber 

bundle  centered  at  4529.2  A)  was  determined  to  be  due 
to  the  scattering  of  the  unshifted  intense  plasma  light. 
Therefore,  it  was  ignored  when  the  ion  emission  from  the 
gap  was  analyzed  to  give  the  electric  field. 


Fig.  2.  (a)  Spectral  profile  of  the  4dg/j  -»  4pj/j  Al++ 
emission  from  the  anode  plasma:  x=0.15  cm  and  1=70 
ns.  The  zero-field  wavelength  is  4529.2  A.  The  initial  A-K 
gap  <fo=0.65  cm  and  B=8.4  kG.  The  inductively  corrected 
dioae  voltage,  monitored  by  a  capacitive  probe  is  shown 
in  the  insert,  (b)  Spectral  profile  from  the  acceleration 
gap  at  x=0.375  cm  and  t=65  ns.  The  line  indicates  the 
trend. 


To  determine  the  electric  field,  wc  calculated  the  the¬ 
oretical  emission  profiles  along  the  magnetic  field  from 
an  ion  in  the  presence  of  perpendicular  electric  Ez  and 
magnetic  B,  fields,  for  the  known  B,  and  a  range  of 
possible  E*.  The  electric  field  was  then  inferred  from 
a  least  squares  fit  of  the  calculated  pattern  to  the  mea¬ 
sured  points  of  the  emission  spectral  distribution.  Since 
this  analysis  assumed  a  single  value  of  the  electric  field,  it 
gives,  for  each  position  x,  the  mean  electric  field  E  over  a 
plane  parallel  to  the  anode. 

Figure  3  shows  examples  of  measured  E(x)  profiles 
for  two  diode  configurations,  one  with  a  planar  cathode 
and  the  other  with  a  vaned  cathode.  In  the  latter,  a 
14-em-long  thin  stainless-steel  vane  was  connected  to  the 
cathode,  parallel  to  the  magnetic  field  lines  projecting  3.5 
mm  into  the  gap  opposite  the  top  of  the  dielectric-anode 
(the  electrons  ExB  drifted  downward).  The  measured 
E(x)  profiles  were  used  to  infer  the  instantaneous  actual 
diode  gap  with  an  uncertainty  of  ±0.075  cm.  As  seen  in 
Fig.  3,  the  field-excluding  anode  plasma  expanded  about 
3  mm  away  from  the  solid  anode  surface  during  the  first 
35  ns  of  the  pulse,  implying  a  plasma  expansion  velocity 
of  5-15  cm/ps.  This  decrease  of  the  diode  gap  results  in 


Fig.  3.  a)  Measured  electric  field  distribution  E(z)  for  a 
planar  diode  configuration  with  d,=0.75.  The  solid  anode 
surface  is  at  x=0  and  the  cathode  surface  is  indicated  by 
the  C.  b)  £(x)  for  a  vaned  cathode  configuration  with 
de= 1  cm;  the  vane  length  /  was  3.5mm.  The  lines  are 
calculated  |4|  E(x)  profiles  using  the  actual  gap. 


an  increase  of  a  factor  of  2-4  in  the  space-charge-limited 
ion  current  density  over  the  value  calculated  using  the 
mechanical  (initially  set)  diode  gap. 

We  have  used  the  actual  diode  gap  to  compare  our 
E(x)  profiles  to  calculations  based  on  the  brillouin-flow 
model  (4).  Such  a  comparison  (see  Fig.  3),  when  done 
for  many  shots,  shows  a  trend:  the  measured  electric  field 
is  higher  than  the  calculated  one  in  the  region  close  to 
the  anode.  This  trend,  which  is  more  pronounced  for 
the  vaned-cathode  configuration,  implies  electron  pres¬ 
ence  beyond  the  theoretical  (4,7)  electron-sheath  region. 

Analysis  of  Diode  Space  Charge 

A  detailed  description  of  the  electron  density  distri¬ 
bution  and  its  implications  inferred  from  our  measured 
E[ x),  will  be  given  elsewhere  (8).  To  summarize,  we  ob¬ 
tain  the  potential  profile  v?(x)  by  <p(x)  =  V  -  f*  E(x)dx, 
and  the  ion  velocity: 

t),(x)  =  (2e(V  -  yj(x))/mi)l/J,  (l) 

where  e  is  the  ion  charge,  is  the  ion  mass,  V  is  the 
instantaneous  diode  voltage,  and  v,(0)=0.  The  ion  cur¬ 
rent  density  J,-,  measured  by  a  Faraday  cup,  is  used  then 
to  give  the  ion  density  n,(x)  =  J,/et>j(x).  We  define  the 
total  electron  number  Ne(x)  =  /*'  ne(x)dx  per  cm 7  be¬ 
tween  the  cathode  plasma  boundary  xc  and  a  position  x 
in  the  gap.  From  Poisson’s  equation,  it  is  given  by 

fVe(x)=5.53  10"E(x)  +  W,(x),  (2) 

where  IV,(x)  is  the  total  ion  number  per  cm5  between  xc 
and  x,  and  E(x)  is  expressed  in  MV/cm.  Using  Eq.  2  we 
infer  the  total  electron  population  Nt  per  cm7  in  the  diode 
gap.  We  then  use  the  model  of  ref.  4  to  calculate  (using 
the  actual  gap)  the  theoretical  thickness  a  of  the  electron- 
sheath,  the  theoretical  ion  current  density  Juh,  and  the 
total  electron  number  Nik,  per  cm7  in  the  sheath.  In  ad¬ 
dition,  vie  obtain  (using  Eq.  2)  the  experimental  electron 
number  Ne,  per  cm7  in  a  slab  of  thickness  s  adjacent  to 
the  cathode  plasma. 

For  the  vaned  cathode  configuration  we  found  that 
Nt/NtK,  and  Nt/N„  are  ~2.  Assuming  space-charge- 
limited  flow  one  would  anticipate  from  this  that  the  ion 
current  density  would  exceed  the  calculated  (4)  one.  In¬ 
deed,  for  this  configuration,  we  observed  Ji/Jhk  between 
2  and  6.  The  uncertainty  is  dominated  by  the  uncertainty 
in  the  inferred  d,  which  again  illustrates  the  need  to  know 
the  actual  diode  gap  in  order  to  understand  the  opera¬ 
tion  of  high-power  diodes.  We  note  that  in  this  analysis 
we  assumed  J,  was  uniform  over  the  entire  anode.  The 
enhanced  electron  presence  near  the  anode  for  the  vaned- 
cathode  configuration  presumably  results  from  the  vane 
tip’s  emitting  electrons  which  approach  the  anode  and 
ExB  drift  in  front  of  it  over  most  of  its  area. 


Transverse  Velocity  Measurements 

Figure  4  shows  measured  profiles  of  the  C  +  +  4647  A 
emission  from  the  anode  plasma  and  from  the  diode  gap. 
The  Zeeman  splitting  for  this  line  is  ~0.2  A  and  its  Stark 
shift  under  our  electric  field  is  <0.1  A.  Thus  the  broad¬ 
ened  line  profile  in  the  gap  gives  the  distribution  of  the 
C++  ion  velocities  parallel  to  the  magnetic  field.  In  Fig.  4 
we  plot,  for  example,  a  Gaussian  fit  to  the  data.  However, 
to  estimate  the  uncertainty  resulting  from  the  unknown 
line-shape,  we  have  analyzed  our  data  assuming  Gaussian 
and  Lorentzian  profiles.  Including  the  measurement  er¬ 
ror  we  obtain  that  50%  of  the  C++  ions  have  transverse 
velocities  less  than  (4  ±1.5 )cm/fia.  Using  the  ion  axial  ve¬ 
locity  at  the  measurement  axial  position,  185 em/jis  (from 
Eq.  1),  we  conclude  that  the  divergence  angle  of  50%  of 
the  C++  ions  is  less  than  (1.4°  ±  0.5°)  in  the  direction 
parallel  to  B  at  the  x  position  of  the  measurement. 

Using  the  3602  A  A1++  line  we  found  similarly  that 
50%  of  the  AI++  ions  have  transverse  velocities  less  than 
(3.4  ±  1.0 )cm/ps,  which  corresponds  to  a  divergence  an¬ 
gle  of  less  than  (1.6°  ±  0.5°),  again  at  the  x  position  of 
the  measurement.  We  also  compared  best  fits  of  the  same 
function  (Gaussian  or  Lorentzian)  for  the  A1  +  +  and  the 
C++  data.  This  comparison  suggests  that  the  transverse 
velocities  scale  as  (e/m,)1/7.  However,  transverse  veloci¬ 
ties  not  dependent  on  e/m,  was  also  consistent  with  the 
data.  Since  our  measurements  were  made  at  a  diode  po¬ 
sition  at  which  the  ions  had  reached  1/3  to  1/2  their  full 
energy,  if  they  do  not  gain  additional  transverse  veloc¬ 
ity  before  leaving  the  diode,  their  divergence  angle  at  ex¬ 
traction  would  be  about  2/3  of  the  values  quoted  above. 
These  are  to  be  compared  with  previous  measurements  (6| 
of  the  proton  divergence  angle  outside  the  diode  parallel 
to  B  which  yielded  4-5°. 
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Fig.  4.  The  open  circles  constitute  the  spectral  profile  of 
the  C++  4647.4  A  line  emission  from  the  anode  plasma 
(x  =  1  mm,  t  -  50  ns).  Curve  (a)  and  the  solid  points  give 
the  system  response.  The  emission  from  the  gap  (t  =  5.5 
mm)  averaged  over  t  =  30  -  70  ns  is  shown  by  +  points 
and  curve  (b)  is  a  Gaussian  best  fit. 


Conclusions 

The  results  described  here  from  the  electric  field  mea¬ 
surement,  namely  the  rapid  early  time  gap  closure  and 
the  presence  of  substantial  electron  space  charge  near  the 
anode,  especially  in  diodes  with  cathode  vanes,  provide 
an  explanation  for  the  enhanced  ion  current  density  in 
magnetically  insulated  diodes  operating  at  modest  power 
densities  such  as  ours.  Taken  together  with  the  highly  dia¬ 
magnetic  electron  flow  characteristic  of  high  power  den¬ 
sity  diodes,  as  discussed  in  ref.  9,  our  results  may  also  ex¬ 
plain  the  enhancements  of  40  or  more  discussed  by  Miller 
(10|  in  experiments  on  Proto  I,  Proto  II  and  I’BFA  I. 

Finally,  we  note  that  if  the  diagnostic  tool  described 
here  is  used  to  obtain  information  about  the  charge  flow 
and  the  ion  transverse  velocities  in  diodes  with  higher 
electric  field,  the  Stark  shift  and  the  Doppler  broadening 
(for  a  given  divergence  angle)  will  be  larger.  Moreover, 
with  higher  electric  fields,  there  will  be  a  higher  ion  cur¬ 
rent  density,  and  lower-lying  ionic  energy  levels  will  be  suf¬ 
ficiently  shifted  to  be  useful,  both  of  which  will  contribute 
to  having  more  signal.  Therefore,  the  measurements  de¬ 
scribed  here  should  be  easier  in  diodes  with  higher  electric 
fields  than  the  one  used  here. 
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A  novel  technique  for  the  production  of 
medium  or  high  energy  electron  beams  from  a  high 
current,  aub  megavolt  pulae  power  drive  electron 
beam  ia  described. 

A  high  current  weakly  relativistic  electron 
beam  ia  used  to  excite  high  power  microwave 
radiation  in  a  transverse  magnetic  mode.  The 
microwave  radiation  is  propagated  through  an 
undulating  guide  which  serves  to  couple  the  fast 
electromagnetic  wave  to  the  electrons.  The 
period  of  the  undulation  ia  chosen  such  that  the 
electrons  slip  one  cycle  of  the  rf  wave  every 
period  of  the  undulation.  The  process 
envisioned,  which  produce*  a  continuous 
acceleration  of  the  electron*,  i*  similar  to 
chat  proposed  in  the  Inverse  Free  Electron  Laser 
accelerator  Cl].  Example*  are  presented  showing 
the  trad*  off  between  the  average  acceleration 
electric  field  gradient  and  the  average  beam 
current. 

Introduction 

V*  examine  in  this  paper  a  novel  technique 
for  the  acceleration  of  electrons  to  medium  or 
high  energies.  The  proposed  technique  uses  a 
high  current  weakly  relativistic  electron  beam 
to  generate  a  high  power  microwave  signal  in  a 
transverse  magnetic  mode  of  the  guide.  These 
processes  have  been  developed  over  the  last 
decade  [2.  3]  and  have  been  shown  to  be  capable 
of  producing  gigawatt  power  pulse*  *t  frequen¬ 
cies  in  excess  of  10  GHs.  The  efficiency  of  the 
microwave  generation  approaches  20J.  The  radia¬ 
tion  is  then  used  to  accelerate  a  lower  current 
secondary  beam  of  electron*  to  high  energy.  The 
acceleration  occur*  ia  a  smooth  bore  undulating 
waveguide  with  the  electrons  slipping  on*  cycle 
of  the  rf  wav*  every  period  of  the  undulation. 
By  us*  of  the  TM02  mode  (or  a  higher  order  node) 
w*  can  us*  the  phase  slip  to  advantage  to  chat, 
as  the  electrons  sample  the  different  radial 
position*  of  the  guide,  the  axial  electric  field 
is  always  in  the  teas*  leading  to  acceleration. 
It  is  worth  noting  cartaia  feature*  of  this 
accelerator  concept.  Puls*  power  sourest  are 


used  to  generate  short  pulses  (~  50-100  nsec.) 
of  high  power  microwave  radiation.  In  order  to 
couple  these  effectively  to  a  secondary  beam  for 
acceleration  it  is  necessary  to  use  a  coupling 
structure  with  a  moderately  high  group  velocity. 
In  the  proposed  accelerator  the  velocity  ia  typ¬ 
ically  of  order  0.4c.  a  value  substantially 
larger  than  that  encountered  in  conventional 
r.f.  linacs.  Other  interesting  and  desirable 
features  include: 

i.  The  us*  of  overmoded  guide  so  that  higher 
frequencies  can  be  used  relatively  easily. 

ii.  The  electric  field  peaks  on  axis  and  not  at 
the  boundaries  of  the  structure.  This 
reduces  the  possibility  of  breakdown  at 
guide  boundaries. 

iii.  The  absence  of  transverse  cavities  reduces 
the  possibility  of  besm  break-up  as  a 
result  of  the  excitation  of  unwanted  cavity 
modes.  This  latter  feature  should  permit 
the  us*  of  larger  beam  currents  without 
serious  instability. 

iv.  Wake  field  effects  should  not  be  any  worse 
chan  those  encountered  in  conventional  (low 

wave  structures  and  may  be  better  at  a 
fixed  frequency  since  we  operate  in  a  high 

order  guide  mod*.  Wall  proximity  effects 
should  therefore  be  reduced. 

We  now  describe  the  accelerator  concept  in 
detail.  Low  average  current  beams  may  be 
accelerated  to  very  high  energies  at  gradients 
in  excsss  of  100  MV/m  or,  by  allowing  a  reduc¬ 
tion  in  the  achievable  gradient,  w*  can 
accelerate  beams  with  average  currents  of  order 
of  a  few  hundred  amperes  at  field  gradients  of 
order  10  MV/n.  The  former  accelerator  confi¬ 
guration  is  of  interest  for  high  energy  physics 
applications  while  the  lover  gradient,  high 
current,  device  may  be  used  for  a  variety  of 
applications  including  the  driving  of  free  elec¬ 
tron  lasers. 

Are.Ufitaf  Caactat 

W*  consider  an  electromagnetic  wav*  pro¬ 
pagating  through  an  undulating  guide.  The  undu¬ 
lation  period  is  P  and  its  amplitude  r  .  The 


axis  of  Che  guide  varies  sinusoidally  as  a  func¬ 
tion  of  posicion  according  to 


r(z)  =  r  sin(2ss/P)  (1) 
o 

The  electron  beam  propagates  through  the  guide 
parallel  to  the  z  axis  with  a  velocity  v  s  c. 
For  simplicity  we  cake  r  /P  <<  1.  The  force 
felt  by  an  electron  due  to  the  rf  wave  is 


F(z)  =  e!  J  Ck  r  sin(2tB/P2coe(4it-kz)  (2) 

0  O  C  0 


and  averages  to  zero  unless  the  synchronism  con¬ 
dition 


^e  *  (k  ♦  2nK)c 


(3) 


applies,  where  k  and  K  are  the  wavenumbers  of 
the  wave  and  the  ripple  respectively,  u  the  wave 
frequency,  the  normalized  electron  velocity, 
c  the  speed  of  light  and  n  an  integer,  which  we 
shall  assume  to  be  unity  in  what  follows.  Using 
the  synchronism  condition  Eq.  2  becomes 


F(z)  *  eE  J  [k  r  sin(2«*/P)  ]cos(4iit/P)  (4) 
0  0  CO 

At  z  *  0.  P/2,  and  P  the  electrons  feel  the  full 
axial  field  Eq  of  the  TM  wave.  Mid  way  between 
these  points,  the  electrons  are  located  close  to 
the  second  radial  extremum  of  the  axial  electric 
field  of  the  TM  mode.  The  field  gradient  is 
shown  as  a  function  of  axial  position  in  Fig.  1. 


Relative  electric  field  felt  by  an  electron  as  a 
function  of  position  in  the  undulator.  The 
electric  field  it  normalised  to  that  on  axis  in 
the  TM31  mode  at  the  same  power  level  and  for  a 
fixed  ratio  of  the  wave  to  the  guide  cut-off 
frequency.  The  solid  line  represents  the  TM01 
mode  with  r/a  *  0.8  and  the  dashed  line  the  TM02 
mode  with  r/a  *  0.6. 


Figure  2  shows  the  field  gradient  averaged  over 
a  period  of  the  undulator  for  the  TK02  mode  at 
various  power  levels.  The  average 


Acceleration  field  averaged  over  the  undulator 
period  for  a  TM02  mode  at  35  GHz  as  a  function 
of  rQ/a  at  power  levels  of  1  CM,  500  MU  and  200 
MU  respectively. 

field  gradient  achievable  in  the  TM02  mode  is 
higher  than  that  achievable  in  the  fundamental 
mode  since  an  electron  must  lose  energy  in  pro¬ 
pagation  through  half  the  undulator  period.  The 
axial  electric  field  in  the  TM  modes  increases 
rapidly  as  the  wave  frequency  decreases  towards 
the  cut-off  frequency  for  the  mode.  A  practical 
limit  it  set  on  the  operating  frequency  by  the 
reduction  in  the  wave  group  velocity  on  the  one 
hand  and  the  increase  in  the  guide  losses  on  the 
other.  At  a  fixed  value  of  the  wave  phase  velo¬ 
city  the  achievable  field  at  fixed  power  through 
the  guide  varies  proportionally  with  the  wave 
frequency . 


Acceleration  in  a  Standing  Haxa 

Ue  consider  the  case  of  an  undulating 
guide,  which  satisfies  the  criterion  of  Eq.  3 
for  wave-electron  synchronism,  at  part  of  a  cav¬ 
ity  of  length  L  >>  P.  The  synchronism  condition 
is  identical  to  that  given  earlier  with  now,  the 
wave  number  k  s  tr  p/L  determined  by  the  resonant 
frequency  of  the  cavity  and  the  cut-off  fre¬ 
quency  of  the  guide.  Mote  that  the  synchronism 
condition  may  be  expressed  in  the  alternate 
form. 


P  *  K 


c  K  *  K 


(5) 


where  Ag,  it  the  guide  wavelength  and  the 
wave  phase  velocity.  This  form  is  instructive 
since  it  shows  that  no  special  tuning  of  the 
guide  undulation  period  is  needed  as  the  elec¬ 
tron  energy  increases.  At  synchronism,  the 

electric  field  felt  by  an  electron  in  traversing 
a  period  of  the  undulator  in  the  cavity  is  given 

bT 


ECi)  *  E  J  (  k  r  sin(2?)  JcosC®4-) 
ax  o  c  o  r  L 

coat (k  ♦  2K)z]  (6) 


This  result  ia  plotted  out  in  Fig.  3. 


Fig.  3 

Local  acceleration  field  normalized  to  the  peak 
electric  field  on  exit.  The  cavity  ia  excited 
by  a  TM02  mode  at  33  GHz.  The  cavity  length  ia 
la  (p  *  78)  and  ro/a  *  0.6.  The  guide 
wavelength  haa  been  choaen  to  fce  equal  to  the 
unduletor  period. 

In  Fig.  4  we  ahow  the  acceleration  field  achiev¬ 
able.  averaged  over  the  unduletor  period.  The 
field  gradient  ia  ah  own  aa  a  function  of  the 
non-dimensional  curvature  (r Q/t)  of  the  guide. 
The  axial  field  atrength  vaa  determined  by 
equating  the  power  loaaea  in  the  walla  (aaauming 
copper  walla  are  room  temperature)  to  the  power 
in  from  the  aource.  Input  power  levela  of  100 
and  300  MW  ware  aaaumed. 


Average  acceleration  field  calculated  for  100 
and  300  MW  at  33  GHz. 

Rote  that  for  the  1  a  cavity  length  aaaumed.  the 
wave  energy  tranait  time.  L/v^ ,  ia  only  10  naec. 
The  average  field  gradianta  achieved  in  thia 
caae  are  comparable  to  thoae  abown  in  Fig.  2  but 


with  power  levels  about  one  order  magnitude 
lower.  An  intereating  variant  ~n  the  cavity 
which  combinea  the  advantages  of  the  cavity  mode 
of  operation  yet  retama  the  desirable  features 
of  the  travelling  wave  accelerator  is  the  trav¬ 
elling  wave  resonator  proposed  in  the  fifties 
[4].  In  this  configuration,  a  single  wave  cir¬ 
culates  around  a  racetrack  with  the  reaonance 
condition  being  act  by  the  cavity  length  being 
an  integral  number  of  guide  vavelengths.  In 
this  caae  the  ratio  of  the  acceleration  field  to 
the  axial  field  ia  the  tame  aa  in  the  travelling 
wave  case. 

Accelerator  Parameter. 

We  present  in  thia  section  typical 
accelerators  paraawtera  for  two  applications, 
namely  a  high  energy  electron  accelerator  and 
secondly,  a  high  average  current  device  at 
medium  energy.  Table  1  showa  the  operating 
parameters  for  a  high  energy  accelerator 
operated  in  the  travelling  wave  TM02  mode  at  33 
GHz.  and  in  the  cavity  mode  at  12.5  GHz.  The 
wave  to  cut-off  frequency  is  1.06  corresponding 
to  a  wave  group  velocity  of  0.33c.  The  travel¬ 
ling  wave  module  length  is  6  m  and  the  cavity 
module  length  1  m.  We  assume  that  the  cavity  is 
pumped  for  one  quarter  of  the  filling  time  2Q/u 
at  the  specified  power  level  and  it  ia  assumed 
that  the  beam  drains  the  stored  energy  com¬ 
pletely.  The  average  beam  current  ia  computed 
on  the  baais  of  a  repetition  rate  of  1  kHz  and 
the  peak  current  on  the  basis  of  allowing  a  1.5Z 
spread  in  the  output  beam  energy. 

Table  I. 

™02  ™02 

Travelling  Wave  Cavity  Mode 


Frequency  (GHz) 

35 

12.5 

Power  (MW) 

1000 

500 

Field  on  Axis  (MV/m) 

310 

450 

Acceleration  Field 

(MV/m) 

100 

100 

Unduletor  Period  (cm) 

2.8 

7.2 

Wall  Modulation 

0.6 

0.6 

Humber  of  Particles/ 

4x1 0H 

3xl012 

Bunch 

Average  Current  (pA) 
Synchrotron  Badiation 

64 

500 

Energy  Loss  (MeV/m 

at  100  CeV 

0.8 

2 

Synchrotron  radiation  losses  have  been 
estimated  for  this  device  and  are  tabulated. 
The  losses  arise  from  the  radial  forces.  Since 
che  primary  acceleration  field  is  axial  the 
radiation  losses  are  minimal  and  are  unimportant 
until  energies  of  order  1  TeV  are  achieved.  The 
losses  scale  as  the  axial  field  squared. 

In  Table  II  ve  show  a  case  representing  the 
generation  of  an  high  average  current  modulated 
beam  suitable  for  use.  for  example,  in  an  FEL 
experiment.  Once  again,  ve  assume  a  module 
length  of  1  m.  In  this  case,  the  output  beam 
consists  of  a  train  of  closely  spaced  bunches 
each  of  which  is  of  very  short  duration  and  are 
separated  from  each  other  by  the  wave  transit 
time  through  tha  cavity  (~  10  nsec.).  In  the 

example  given  the  beam  has  an  average  current  of 
100  A.  a  peak  current  of  1.8  kA.  and  an  energy 

spread  of  about  12. 


Table  II. 


Frequency  (GHz) 

5.7 

Power  (GH) 

1.0 

Field  on  Axis  (H7/m) 

50 

Acceleration  Field  (MV/m) 

10 

Undulator  Period  (cm) 

17 

Wall  Modulation 

0.3 

Number  of  Particles/Bunch 

1.5 

x  10 

Max  Bunch  Rate  (sac  *) 

4  x 

109 

Peak  Current  (kA) 

1.8 

Average  Current  (A) 

100 

Concluding  CfflMSBtl 


Propagation  in  the  TM11  mode  through  a  60  cm 
long  undulator  guide,  with  rg/a  -  O.S.  as  a 
function  of  frequency. 

>>  appear  to  be  substantiated.  He  also  have 
results  shoving  propagation  in  the  TMQ2  mode 
vith  no  noticeable  mode  conversion.  This  case 
is  more  difficult  to  test  experimentally  as  ve 
still  excite  some  TM01  mode  simultaneously  vith 
the  TM02  mode.  In  any  event  ve  find,  even  vith 
quite  severe  curvature  of  the  guide,  that  the 
uniform  guide  node  is  essentially  preserved. 

He  conclude  that  the  proposed  accelerator 
configurations  show  promise  for  use  of  giant 
power  microwave  pulse*  in  electron  accelerators. 
The  final  success  of  this  proposed  accelerator 
probably  depends  on. the  achievable  phase  stabil¬ 
ity  of  the  source*. 

AclmwXnilgaxpti 


We  conclude  this  paper  vith  some  comaants 
regarding  possible  limitations  in  tbs  proposed 
accelerator.  One  concern  is  that  the  wavs  pro¬ 
pagation  will  be  very  different  to  that  assumed. 
It  is  clear  that  ve  have  a  periodic  structure 
and  that  this  will  result  in  a  modification  of 
the  expected  uniform  guide  modes.  To  examine 
this  ve  have  fabricated  a  60  cm  long  test  sec¬ 
tion  of  2.3  cm  radius  guide  vith  r  /a  *  0.3  and 

o 

r  1  16  cm.  Preliminary  tost  results  are  given 
in  Fig.  3  which  shows  propagation  in  the  TM11 
mode  from  4  to  14  GHt . 

A  detailed  analysis  of  the  network  analyser  out¬ 
put  yields  a  dispersion  relationship  similar  to 
that  obtained  in  a  section  of  uniform  pipe. 
There  exist*  an  upward  shift  in  the  frequency  of 
up  to  300  MBs  which  we  attribute  to  distortion 
in  the  guide  cross  section  during  the  bending 
process.  Our  qualitative  expectations  that  the 
dispersion  would  not  change  significantly  if  P 
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Abstract 

Light  ion  inertial  confinement  fusion 
research  at  Cornell  includes  two  ion  diode 
development  experiments  which  utilize  the  O.rtTJ 
(1.8MV,  450kA),  40ns  LION  accelerator.  One 
uses  a  Bp  extraction  geometry  magnetically 
insulated  diode  and  has  tw*o  goals:  to  study  ion 
beam  generation  using  an  exploding  metal  film 
anode  plasma  source,  and  to  study  ion  bean 
propagation  and  focussing.  The  second  experiment 
utilizes  the  ^  "barrel"  geometry  and  its  goals 
are:  to  determine  if  this  symmetric  diode  can  be 
run  on  a  "single  sided"  generator,  and,  if  so,  to 
investigate  the  effects  on  diode  operation  of 
pulse  risetime,  the  details  of  the  magnetic  field 
geometry,  and  other  mechanisms  that  affect 
electron  leakage  to  the  anode.  Results  from 
these  experiments  are  described. 

Introduction 

In  this  paper,  report  research  progress 
at  Cornell  University  since  the  last  beams  con¬ 
ference  r 1 1  on  the  generation  of  intense  ion 
beams  for  inertial  confinement  fusion  (ICF).  Our 
effort  involves  two  ion  diode  development 
experiments  using  the  0.8TW  (1.8MV,  450  kA) ,  40 
ns  LION  accelerator,  a  module  of  the  FBFA  I 
machine.  Che  experiment  utilizes  an  extraction 
geometry  magnetically  Insulated  diode  (MID)  and 
has  two  goals,  the  development  of  an  exploding 
metal  film  active  anode  plasm  source  (EMFAAPS), 
and  the  study  of  the  propagation  and  focussing  of 
an  intense  ion  beam  from  an  extraction  diode. 

The  E7CAAPS  is  driven  by  the  first  few  nano¬ 
seconds  of  the  LION  power  pulse  via  a  plasma 
opening  Mitch.  Using  aluwinun  foils,  about 
500A/o«2  of  ions  heavier  than  protons  (and 
200a /cm2  of  protons)  have  been  produced  at  about 
IMeV.  Propagation  and  focussing  experiments  have 
so  far  been  performed  only  with  an  epoxy-filled 
"surface  flashover"  anode  in  place  of  EMFAAPS. 
Preliminary  results  show  a  radial  aiming  error 
that  may  N»  due  to  a  small  electric  field  (about 
lOkV/cm)  in  the  strong  transecve  magnetic  field 


region  near  the  diode  Insulation  coils,  or  to  the 
virtual  cathode  being  tilted  with  respect  to  the 
anode  surface. 

The  second  experiment  utilizes  the  applied 
Bz  barrel  MID  configuration  of  Johnson  et  al. 

(2) .  The  first  goal  of  this  experiment  was  to 
determine  if  this  cylindrically  symmetric  diode, 
which  is  also  symmetric  about  its  midplane,  can 
be  efficiently  operated  using  a  power  source 
(LION)  which  is  not  inherently  symmetric  about 
the  diode  midplane.  Recent  results  in  which 
reasonably  uniform  ion  beam  have  been  produced 
suggest  that  this  can  be  done.  Wb  are  now 
characterizing  the  operation  of  the  diode  and  the 
ion  beam  it  produces  under  "standard  conditions". 
Wb  will  then  investigate  the  effects  on  diode 
operation  of  pulse  rise  time  and  mechanisms  which 
affect  the  electron  leakage  to  the  anode  (e.g.  - 
an  electron  limiter  (3),  azimuthal  asymmetries, 
and  the  overall  magnetic  field  geometry). 

Extractor  Diode  Experiments 

EMFAAPS 

The  development  of  an  adequate  anode  plasma 
ion  source  for  light  ion  ICF  is  a  high  priority 
problem  at  the  Ignition  experiment  level.  Wb  are 
investigating  an  exploding  metal  film  active 
anode  plasma  source  (EMFAAPS) .  Wb  have  found 
that  a  metal  film  deposited  an  a  plastic  surface 
can  be  made  to  explode  with  reasonable 
uniformity,  as  determined  fra*  the  visible  light 
image  using  a  streak  camera,  when  driven  by  the 
first  few  nanoseconds  of  the  LION  accelerator 
pcwer  pulse  f  4 ] .  This  technique  has  now  been 
tested  using  an  evaporated  metal  film  (500  4  of 
alumlnus)  on  an  epoxy  snode  mounted  in  the  diode 
shear:  in  Fig.  1.  The  film  is  vaporized  and 
undergoes  electrical  breakdown  by  passing  the 
generator  current  through  the  film  via  the  plasma 
switch  shown  inside  the  anode  in  Fig.  1.  High 
voltage  appears  across  the  anode-cathode  gap  as 
the  plans  switch  ixpedance  rises,  allowing  a 
high  current  ion  beat  to  be  extracted.  The 


plasma  switch  geometry  used  here  does  not  open; 
instead  the  switch  current  remains  nearly 
constant  at  80  -  150A  (depending  on  the  precise 
plasna  switch  conditions)  durinq  the  high  voltage 
phase  of  the  pulse,  (see  dashed  trace  in  Fig  2b). 

The  earliest  tests  with  EMFAAPS  produced 
very  low  impedance  diode  operation  and  a  great 
deal  of  electron  damage  to  the  inner  portions  of 
the  anode,  even  though  B/B*  >1.8,  where  B*  is  the 
critical  field  for  magnetic  insulation  f5]. 

Field  mapping  showed  flux  surfaces  crossing  the 
epoxy  anode  surface  at  the  location  of  the 
damage.  The  epoxy  thickness  was  modified  to 
eliminate  the  intersections  of  flux  surfaces  and 
the  epoxy  surface  resulting  in  the  diode 
electrode  and  flux  surface  oonf igu rat  ions  shown 
in  Fig.  3.  The  result  was  decreased  electron 
damage,  and  higher  inpedance  diode  operation, 
as  shown  in  Fig.  2.  In  this  shot,  peak  voltage 
was  about  1  MV  and  peak  total  current  was  about 
400  kA.  Of  that  total,  about  !50kA  is  switch 
current,  150  kA  (peak)  is  ion  current  and  the 
remainder  is  electron  leakage  current. 

Also  shown  in  Fig. 2  are  three  Faraday  cup 
traces  frcm  cups  located  9  and  14  an  fran  the 
anode.  Two  current  peaks  are  visible.  The  first 
is  consistent  with  prompt  tum-on  of  the  anode  if 
the  time-of-f light  of  protons  is  taken  into 
account.  The  second  is  consistent  with  pronpt 
anode  tum-on  and  time-of-f  light  for  A1+.  The 
proton  contamination  is  to  be  expected  since  no 
effort  was  made  to  avoid  it.  Thus,  although  we 
have  not  confirmed  that  the  heavier  species  is 
A1+,  we  believe  that  we  are  seeing  a  peak  current 
density  of  400-700  A/cm 2  of  that  species  a 
distance  of  9  am  fran  the  anode.  Since  the  area 
of  the  anode  is  about  160  cm2,  we  estimate  a 
total  peak  A1+  current  of  (90i30)  kA.  Because  of 
the  time  of  flight  spreading  of  the  beam  in 
propagating  to  the  Faraday  cups,  this  may  be  an 
underestimate.  The  proton  peak  current  density 
of  250-400  A/cm2  similarly  extrapolates  to  a 
total  peak  proton  current  of  (50*10)  kA. 
tkxever,  notice  that  the  proton  peak  Is  only 
about  15  ns  wide,  suggesting  that  the  proton 
component  in  the  anode  plaena  is  depleted  in  the 
first  few  ns  of  the  power  pulse  using  EMFAAPS. 

After  each  shot,  the  aluminum  film  is 
conpletely  gone  and,  in  the  less  damaging  pulses, 
electron  damage  is  aonflned  to  a  snail  region  of 
the  epoxy  substrate  near  the  Inside  radius  of  the 
anode. 


Ion  Beam  Propagation  and  Focussing  Experiments 

The  diode  shown  in  Fig.  1  has  also  been  used 
with  a  more  conventional  epoxy-f illed-groove 
aluminum  anode  instead  of  the  EMFAAPS  to  extract 
an  ion  beam  and  study  its  propagation  and 
focussing.  Measurements  of  the  mean  direction  of 
motion  and  the  intrinsic  divergence  of  a  portion 
of  the  ion  beam  have  been  made  by  placing  a 
pinhole  aperture  9  an  from  the  anode  (just 
behind  the  cathode  flange  in  the  field  froe  drift 
space,  see  Fig.  1).  The  selected  beamlet  is  then 
incident  on  an  array  of  charge  collectors,  the 
geometry  of  which  is  shown  in  Fig.  1.  Because  of 
the  short  (50  ns)  diode  power  pulse  and  the 
nearly  20  ns  tum-on  delay  with  epoxy-f  illed- 
groove  anodes,  time-of-f light  separation  allows 
the  unfolding  of  direction-of-motion  information 
for  the  various  ion  species  in  the  beam,  and 
divergence  information  for  the  protons. 

Figure  4  shows  data  from  a  shot  with  the 
aperture  located  near  the  outer  edged  of  the  beam 
envelope.  Ions  began  to  appear  at  about  30  ns 
and  the  peak  of  the  Ion  power  pulse  was  at  about 
53  ns  on  this  shot.  At  that  time  the  proton 
divergence  is  about  3  degrees  FVMM.  The  peaks 
for  ions  heavier  than  protons,  assumed  to  C+  and 
C++,  overlap  for  much  of  the  pulse  so  determining 
the  divergence  for  these  lone  is  not  possible. 
During  the  pulse,  the  ratio  of  C++  to  C+  seems  to 
increase  although  again  because  of  the  overlap  it 
is  difficult  to  quantify  the  ratio.  The  proton 
and  the  C++  peaks  move  rapidly  Inward  vith  time. 
Figure  5  shows  a  oonparison  of  the  observed 
deflection  of  the  ions  with  a  calculated 
deflection  that  accounts  for  magnetic  bending  but 
assume  there  is  no  space  charge  electric  field 
(l.e.  -  the  beam  is  fully  charge  neutralized). 

The  calculation  assumes  that  lone  leave  the  anode 
normal  to  the  surface  with  an  energy  determined 
by  the  instantaneous  value  of  the  voltage,  that 
the  beai  current  Is  neutralized  once  it  passes 
through  the  entrance  of  the  flux  conserving 
cathode  extraction  flange,  and  that  canonical 
angular  momenta*!  is  conserved.  The  actual 
deflection  of  the  proton  oowponent  of  the  bean  Is 
consistently  2-4  degrees  less  than  the  calculated 
values.  Figure  5  also  shows  the  measured  and 
calculated  deflection  for  C+  tons,  *rtrich  are  not 
as  affsctsd  by  the  magnetic  fields.  If  we 
estimate  the  space  charge  electric  field  required 
to  deflect  the  ions  by  l.S*  assuming  it  is 


present  uniformly  whereever  there  is  strong 
maqnetic  field,  we  find  that  a  field  of  about  10 
kV/cm  is  enough.  This  corresponds  to  >  99.9% 
charge  neutrality  of  the  ion  beam.  Although  the 
ion  mean  deflection  angle  was  also  obtained, 
and  the  results  regarding  a  space  charge  electric 
field  obtained  with  them  is  consistent  with  the 
above  discussion,  we  do  not  include  them  in  Fig. 

5  in  order  to  avoid  overly  conplicating  it. 

Another  possible  source  of  the  observed 
outward  deflection  of  the  ions  is  a  virtual 
cathode  which  is  not  parallel  to  the  anode 
surface.  This  might  result  from  the  inherent  1/r 
dependence  of  the  insulating  magnetic  field.  If 
the  virtual  catliode  is  tilted  and  tl>e  ion  current 
density  is  space  charge  limited,  then  the  ion 
current  density  should  give  an  indication  of  the 
variation  of  the  effective  gap  spacing  d  through 
the  l/d?  current  density  dependence.  Figure  6 
shows  signals  from  three  Faraday  cups  radially 
separated  by  0.8  cm,  placed  inmediately  behind 
r»ie  cathode.  Notice  that  the  highest  current 
)?nsity  is  seen  by  the  outer  cup,  the  middle  one 
is  next,  and  the  inner  one  is  lowest.  Figure  6 
also  shews  the  values  of  d  obtained  using  the 
space  charge  limited  current  density  formula, 
with  the  instantaneous  voltage  taken  from  the 
upper  trace  in  the  figure.  The  difference  in  d 
implied  by  the  detectors  taken  in  pairs  divided 
by  the  spacing  between  them  gives  an  angle  (in 
radians)  which  may  be  Interpreted  as  the  angle  of 
the  virtual  cathode  with  respect  to  the  anode 
surface.  This  angle  is  shown  in  degrees  in  the 
bottom  graph  in  Fig.  6.  Since  the  curves  are 
constructed  from  the  difference  between  two  noisy 
signals,  the  details  of  the  two  curves  are 
probably  not  significant.  However,  the  general 
trend  of  this  angle  would  imply  an  outward 
deflection  of  the  ions  early  in  time  which  would 
gradually  decrease  to  no  deflection  at  the  end  of 
the  pulse,  as  seen  for  the  carbon  ions  (Fig.  5) . 

B^Dlode  experiments 

The  diode  used  here  (see  Fig  7)  is  a 
replica  of  the  Proto  I  ion  diode  used  by  D.J. 
Johnson  et  al.  |2|  at  Sandia  but  modified  as 
needed  to  be  used  on  the  LION  accelerator.  The 
cathode  is  made  of  aluminum  and  is  a  flux 
excluder.  In  contrast  to  nultimodule  pulsed 
power  generators  such  as  Proto  I  and  the  PBFA  I 


and  II  halves  of  the  Bz  diode  symmetrically,  here 
we  mist  provide  paths  to  the  two  sides  of  the 
diode  so  as  to  balance  the  current  flow  from 
LION.  The  geometry  shorn  in  Fig.  7  is  designed 
to  have  equal  inductances  for  the  two  current 
paths.  Experiments  presently  in  progress  have 
determined  that  it  is  possible  to  balance  the 
feeds  and  generate  a  uniform  ion  beam  at  a  diode 
voltage  up  to  about  1  MV.  Figure  8  shows  a 
sample  set  of  data  from  a  shot  in  which  a  damage 
target  showed  that  a  reasonably  uniform  ion  beam 
was  produced.  On  this  shot  the  peak  inner  and 
outer  currents,  I,,  and  I2,  were  llOkA  and  120kA 
respectively,  with  an  incident  current  of  260kA. 
Leakage  current  at  the  post-hole  convolute  was 
measured  to  be  <  5%  of  the  incident  current. 

At  the  present  time  we  are  developing 
techniques  to  characterize  diode  operation  and 
ion  beam  quality,  and  we  will  shortly  undertake 
the  following: 

(1) .  Investigate  the  effect  on  the  turn-on  delay 
of  an  epoxy-f illed-groove  flashover  anode  as  the 
pusle  rise  time  is  varied  with  a  plasma  opening 
switch. 

(2) .  With  a  gas  cell  installed  in  the  volume 
inside  the  cathode,  study  the  current 
neutralization  phenomenon  in  detail. 

(3) .  Determine  the  effect  of  anode  geometry  and 
magnetic  diffusion  into  the  anode  on  the  ion 
focussing. 
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Ahetract 

A  two  stage  proton  induction  linac  has  been 
fabricated  and  the  beam  characteristics  meaa- 
ured.  The  injected  beam  has  an  energy  of  1  MeV 
with  an  ion  current  of  up  to  6.0  kA.  The  post 
acceleration  gap  increments  the  energy  by  350 
keV. 

In  thia  paper,  we  shall  review  the  genera¬ 
tion  and  transport  of  the  proton  beam  through 
the  accelerator.  Time  resolved  ion  profiles 
will  be  presented.  The  dynamics  of  the  ion  beam 
will  be  reviewed  and  examined.  Time  integrated 
proton  flux  and  divergence  sMaaurements  will 
alto  be  presented  as  a  function  of  position 
throughout  the  accelerator. 

We  shall  also  review  the  possibility  of 
using  space  charge  neutralized  ion  beama  in  mul¬ 
tistage  accelerators. 

Introduction 

We  consider  in  thia  paper  the  uae  of  linear 
induction  accelerators  for  the  production  of 
high  current,  multi-magavolt  ion  beams  Cl, 2, 3}. 
For  beam  currents  of  several  kilo-Amperes  the 
self  electric  field  of  the  beam  ia  of  order  50 
KT/a  and  cannot  be  confined  by  modest  (-1  Tesla) 
magnetic  fields.  In  this  experiment  we  present 
data  on  the  design  and  performance  of  a  two 
stage,  6  kA,  proton  induction  linac.  The  ion 
beam,  which  ia  magnetically  confined,  ia  space 


charge  neutralized  by  co-drifting  electrons  ori¬ 
ginating  from  the  electron  cloud  defining  the 
cathode  plane.  Additional  neutralization  is 
provided  by  electrons  emitted  from  the  drift 
tube  wall. 

Pargitigfrg 

The  two  stage  linear  induction  accelerator 
has  an  inductively  driven  first  stage  which 
gives  an  injection  energy  of  1  MeV  and  a  second 
stage  post-acceleration  of  350  kV;  the  half 
width  of  the  pulse  is  60  ns.  Proton  currents  of 
up  to  6  kA  have  been  produced.  The  beam,  which 
is  generated  in  a  magnetically  insulated  diode, 
is  propagated  in  a  solenoidal  magnetic  field  to 
a  cusp  magnetic  insulating  field,  where  the  post 
acceleration  occurs.  The  experimental  confi¬ 
guration  ia  illustrated  in  Fig.  1. 

EC! *  Diagnostics 

Measurements  were  made  of  the  ion  and  net 
current  densities  using  radial  arrays  of  Faraday 
cups.  In  order  to  measure  the  ion  current  den¬ 
sity.  the  electrons  co-drifting  in  the  beam  must 
be  separated  out.  This  is  accomplished  by  the 
detector  illustrated  in  Fig.  2,  The  axial  mag¬ 
netic  field  which  confines  the  beam  is  excluded 
from  the  aperture  in  an  iron  plate.  The  beam 
electrons  follow  the  field  lines  away  from  the 
aperture  t&ile  the  ions  pass  through  the  aper¬ 
ture  to  the  collector.  As  a  check  on  the  per- 


Fig.  1.  Experimental  configuration:  (a)  first  accelaration  gap,  (b)  anode,  (c)  proton  emission 
surfacs,  (d)  electron  amission  fins,  (e)  field  coil,  (f)  solenoid,  (g)  magnetic  field  lines,  Cb) 
second  acceleration  gap,  (i)  fins. 
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—collector 


— magnetic 
field  line* 


—iron  plate 

Fig.  2.  Ion  Faraday  collector.  The  magnetic 
field  line#  are  only  shown  entering  the  iron 
plate. 

foraance  of  the  detector.  the  time  integrateu 
signal  was  compared  with  a  carbon  activation 
measurement.  The  Feraday  cup  signal  was  59  T  8  ! 
of  the  carbon  signal.  This  difference  ia  pri¬ 
marily  attributed  to  secondary  electrons  emitted 
from  the  inside  surface  of  the  aperture  and  col¬ 
lected  by  the  cup. 

The  net  current  Faraday  cup  differs  from 
the  ion  collecting  cup  in  that  the  front  aper¬ 
ture  plate  ia  brass  and  does  not  limit  electron 
flow  to  the  cup. 

Carbon  activation  was  used  to  measure  the 
beam  fluence  over  areas  large  compared  to  the 
1.6  am  aperture  of  the  Faraday  cups. 

Pyp-r  -iiMnr  a  1  B»«nl  ta 

Beam  IraaipOEt. 

Carbon  activation  amesurements  of  the  beam 
fluence  are  shown  in  Fig.  3  as  a  function  of 
radius  for  several  axial  positions.  The  beam 
expands  upon  entry  into  the  drift  tuba.  and  is 
in  contact  with  the  drift  tube  well.  The  beam 
then  detaches  from  the  wall  and  moves  inward • 
where  it  propagates  with  little  further  loss. 
The  beam  does  not  follow  single  particle  trajec- 
tories.  which  have  a  calculated  focal  length  of 
35  cm.  At  e  distance  of  46  cm  from  the  diode* 
the  beam,  while  densest  on  axis,  has  cmly 
focussed  to  a  swan  radius  of  3.3  cm  from  an 
average  radius  of  5. 5  cm  at  injection. 

Diodes,  similar  to  that  used  in  this  work, 
typically  produce  ion  beams  with  a  mean  diver¬ 
gence  angle  of  about  3*.  Measurements  in  this 
sxperiment  indicate  that  the  outermost  protons 
expend  radially  at  a  14*  angle  prior  to  entering 
the  drift  tuba  2.5  cm  downstream  of  the  cathode 
plane.  The  subsequent  radially  inward  proton 
motion  indicates  that  the  ion  beam  is  nautral- 


Fig.  3.  Beam  fluence  aa  a  function  of  radius  at 
(a)  3  cm.  (b)  17  cm.  (c)  32  cm.  and  (d)  44  cm 
from  cathode. 


ized  by  electrons  emitted  from  the  drift  tube 
wall.  Approximately  20X  of  the  diode  generated 
beam  is  lost  to  the  support  struts  for  the  inner 
magnetic  field  ooil  and  a  further  5X  to  the  wall 
at  the  entrance  to  the  drift  tube.  Following 
the  beam  detachment  from  the  wall  the  transport 
through  the  first  section  of  drift  tube  is  100X 
efficient. 

A  net  current  Faraday  cup  was  used  to  sam¬ 
ple  the  beam  as  n  function  of  distance  from  the 
diode.  The  traces,  shown  in  Fig.  4,  indicate 


Fig.  4.  Net  currant  density  at  0.6  cm  ia).  3.5 
cm  Cb),  7.5  cm  (c)  and  12.5  cm  (d)  from  cathode. 
The  detector  was  at  5  cm  radius. 


nanoseconds 

Fig.  5.  Ion  current  density  at  1.2  and  3.5  cn 
radial  positions  for  first  (a)  and  (b)  and  sub* 
sequent  shots  (c)  and  (d). 

that  the  net  current  density  decreases  froa  40 
to  5  A/cm2  over  the  first  12.5  ca  froa  the 
diode.  This  is  consistent  with  the  neutraliz* 
ing  electrons  entering  the  beaa  along  aagnetic 
field  lines  at  the  entrance  to  the  solenoid. 

(See  Fig.  1) . 

Additional  confirmation  of  the  importance 
of  the  walls  in  the  space  charge  neutralization 
process  was  provided  in  a  series  of  shots  util¬ 
izing  two  irot  aperture  Faraday  cups  located  at 
radii  of  1.2  and  3.5  ca  and  at  an  axial  location 
30  ca  downstream  of  tha  diode.  Figure  5  shows 
the  output  froa  the  two  Faraday  cups  during  a 

sequence  of  shots.  The  upper  traces,  taken  on 
the  first  shot  in  the  series,  show  peak  ion 

current  densities  of  12  and  76  A/ca2  on  the 
inner  and  outer  cups  respectively.  The  second 
set  of  data,  taken  two  shots  later,  show  ion 

current  densities  at  the  cups  of  80  and  20  A/ca2 
respectively.  The  current  densities  stabilised 
at  the  second  set  of  values  until  the  tube  wall 
was  cleaned,  after  which  the  signals  returned  to 
their  earlier  values.  Tha  gross  difference  in 
beaa  focussing  seen  here  was  due  to  additional 
beaa  space  charge  neutralization  arising  froa 
the  debris  due  to  the  blow  off  of  dry  silver 
copy  paper,  used  in  the  experiment  to  generate 
an  qualitative  image  of  the  time  integrated  beaa 
profile.  The  effects  of  tha  wall  contamination 
did  however  provide  strong  evidence  for  the 
importance  of  electron  emission  froa  the  walls 
on  the  beaa  neutralization  and  dynaaics,  and 
night  be  used  to  advantage  in  the  design  of 
electron  sources  for  space  charge  neutralisa¬ 
tion. 

Data  froa  Faraday  cups  can  be  used  to  fora 
radial  profiles  of  the  current  density.  Pro- 


(c)  (d) 


Fig.  6.  Sadial  ion  current  density  profiles, 
(a)  40  ns,  (b)  60  ns,  (c)  80  ns  (d)  100  ns  after 
initiation  of  diode  voltage. 

files  at  three  different  axial  locations  are 
shown  in  Fig  6.  Note  that  the  beaa  moves  inward 
about  70  ns  after  the  onset  of  the  diode  voltage 
pulse.  The  delay  is  due  in  part  to  the  ion 
transit  tina  through  the  drift  tube  and  also  to 
delays  in  the  onset  of  electron  neutralization 
of  the  beaa  space  charge.  ' 

The  signals  from  the  iron  and  brass  aper¬ 
ture  Faraday  cups  can  be  coabined  to  estimate 
the  beaa  potential  and  tha  degree  of  neutraliza¬ 
tion  in  the  beaa.  The  difference'  between  the 
signals  is  tha  electron  current  density.  Assum¬ 
ing  that  the  ion  and  electron  charge  densities 
are  nearly  equal,  we  estimate  the  velocity  of 
the  co-drifting  electrons  froa 


The  ion  velocity  is  determined  froa  the  diode 
voltage.  Since  the  electrons  are  ositted  froa 
the  cathode  at  ground  potential,  their  velocity 
provides  a  measure  of  the  local  potential  inside 
the  beaa. 


nanoseconds 

Fig.  7.  Beaa  potential  at  3  ca  radius,  27  ca 
froa  cathode.  The  tiae  scale  is  not  referenced 
to  the  beaa  pulse. 


Figure  7  shove  an  example  of  the  potential,  cal¬ 
culated  in  this  sinner,  at  a  3  cm  radius,  27  cm 
from  the  diode.  The  space  potential  varies 
radially  and  on  a  shot  to  shot  basis.  Values  up 
to  a  factor  of  three  larger  than  thsc  shown  in 
the  figure  have  been  recorded.  The  electron 
drift  velocity  calculated  using  Eq.  1  gives  a 
velocity  approximately  equal  to  or  slightly 
larger  chan  the  ion  drift  velocity.  This  result 
is  vichin  the  range  expected  on  the  basis  of  the 
one  dimensional  theory  of  space  charge  neutrali¬ 
zation.  Finally,  ve  use  the  space  potential  to 
check  on  che  assumption  of  equal  ion  and  elec¬ 
tron  densities.  A  characteristic  solution  of 
Poisson's  equation  shows  che  charge  imbalance  to 
be  of  order  of  0.02Z,  and  in  accordance  with  the 
assumptions  made  in  deriving  the  space  poten¬ 
tial. 

Post  Acceleration 

3eam  propagation  has  been  measured  across 
the  second  gap  using  carbon  activation. 

At  currents  above  1.3  kA.  the  protons  in 
the  beam  will  slow  down  upon  entry  into  the  gap. 
increasing  che  local  charge  density.  In  order 
to  minimize  the  divergence  caused  by  space 
charge  in  che  gap  the  beam  current  was  decreased- 
below  this  value.  The  propagation  was  measured 
by  four  targets  just  upstream  of  the  second  gap. 
The  beam  fluence  upstream  could  then  be  calcu¬ 
lated  by  scaling  of  Che  proton  fluence  incident 
on  these  targets;  the  downstream  intensity  was 
measured  over  the  entire  beam  area.  The  carbon 
targets  were  placed  in  an  annular  pattern  so 
that  the  radial  profile  of  the  beam  could  be 
measured  in  addition  to  the  total  transport. 
Beam  fluence  profiles  upstream  and  downstream  of 
the  second  gap  are  shown  in  Pig.  8.  Propagation 
efficiency  was  measured  at  88Z  for  a  1500  A  shot 
sad  100Z  for  an  iOO  A  shot.  The  standard  devia¬ 
tion  error  vs a  51. 

He  have  alto  used  the  radial  profiles  to 
obtain  an  estimate  of  the  beam  divergence  in 
propagation  through  the  second  gap.  He  consider 
chose  protons  incident  at  che  gap  between  radii 
of  3. A  and  4.7  cm  and  calculate  the  fraction  of 
those  protons  which  finish  downstream  of  the  gap 
in  the  radial  interval  between  4.7  sad  6.0  cm. 
The  gap  width  it  3.8  cm.  He  find  36Z  and  31Z  of 
the  protons,  corresponding  to  the  1300  sad  800 
Ampere  beams  reapectively,  make  this  transition. 
For  comparison,  we  note  that  unnautraliaed  and 
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Fig.  8.  Seam  fluence  (a)  upstream  of  the  second 
gap;  (b)  downstream  of  the  second  gap  (1500  A); 
(c)  downstream  of  the  second  gap  (800  A). 

uncoax ined  beams  of  the  same  diameter  and  carry¬ 
ing.  the  same  currents  as  those  reported  above 
would  have  50Z  and  31Z  of  the  protons  make  the 
radial  transition. 

Conclusions 

He  have  developed  and  used  an  array  of 
diagnostic  techniques  co  study  beam  transport 
and  neutralization  along  a  uniform,  strong  axial 
magnetic  field.  The  same  techniques  have  been 
used  to  study  beam  transport  through  a  350  kV 
post  acceleration  gap.  The  beam  dynamics 
observed  approximate  those  expected  on  che  basis 
of  tingle  particle  orbits  in  the  field  geometry 
employed.  An  additional  beam  divergence  it 
identified  vith  che  delay  associated  with  space 
charge  neutralization  of  the  beam  channel.  Very 
efficient  beam  transport  has  been  observed 
through  the  post  acceleration  diode.  The  beam 
divergence  in  that  section  of  channel  will 
require  additional  focussing  in  the  gap,  prob¬ 
ably  through  field  contouring. 
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Abstract 

The  propagation  of  intense  ion  beams  in  an 
ionized  medium  and  z-pinch  current  channels  is 
examined.  The  analysis  of  beam  dynamics  includes 
both  the  induced  return  current  and  the  self- 
consistent  nonlinear  response  of  the  plasma 
ions.  Stable  propagation  in  channels  is  possible 
for  regimes  of  interest  to  inertial  fusion. 

Introduction 

An  understanding  of  intense  ion  beam  propa¬ 
gation  is  particularly  inportant  for  inertial 
fusion  reactor  studies,  and  for  injection  of  ion 
beams  into  magnetic  traps  for  magnetic  fusion. 

It  may  also  have  some  relevance  to  the  propaga¬ 
tion  of  plasmoids  in  the  solar  atmosphere  result¬ 
ing  from  solar  flares,  to  cometary  tails,  and 
possibly  to  such  astrophysical  phenomena  as 
jets.  This  paper  addresses  the  propagation  of 
intense  ion  beams  in  (i)  an  ambient  ionized 
medium  and  (il)  in  preformed  z-pinched  current 
channels  on  a  long  enough  time  scale  to  include 
the  motion  of  the  background  plasma  Ions.  A 
brief  discussion  of  the  expected  major  instabil¬ 
ities  is  also  undertaken. 

Early  experiments  [1]  confirmed  that  intense 

ion  beams  propagated  freely  in  vacuun,  provided 

electron  emitting  surfaces  were  available  nearby 

for  the  almost  complete  space  charge  and  current 

neutralization  that  is  necessary.  An  elementary 

one-dimensional  calculation  (2]  accounts  for 

electron  acceleration  frcm  such  a  surface  ty  a 

potential  *Q  =*  2m0v^/e  that  develops  across  a 

narrow  sheath;  v  is  the  beam  velocity,  and  e 
o 

and  m  are  the  electron  charge  and  mass,  respec¬ 
tively.  C.n.S.  electrostatic  units  will  be  used 
throughout.  The  accelerated  electrons  eventually 
thermallze  through  collective  interactions  13,4) 
such  that  they  acquire  a  thermal  velocity  vq 
over  and  above  their  mean  velocity  v  with  a  mean 
bean  potential  eT  »  1/2  m@v£. 

When  a  charge  neutralized  non-relativlstic 
Intense  ion  beam  is  injected  across  a  static 
magnetic  field  Bz  in  vacuun,  it  may  propagate 
only  if  a  polarizing  electric  field  *  vqB/c 
can  be  maintained  across  the  bean  thickness  by 


virtue  of  the  y  displacement  of  the  acconpanying 
electrons.  This  is  possible  if 


4nn.  m.  AS?  4C  1 

b  b  o 

(la) 

(4wn. m.  c^/B^)  (a/o)  »  1 
b  b 

(lb) 

are  satisfied  (5,6);  n^n^is  the  mass  density  of 
the  beam,  a  is  the  half  thickness  and 
p  =  m^v^c/eB^  is  the  beam  ion  Larmor  radius,  the 
subscript  b  refers  to  bean  quantities  throughout. 
When  E  is  shorted  out  in  a  time  less  than  the 

y 

pulse  duration  t  either  by  the  magnetic  geometry 

or  background  plasma,  then  propagation  ceases. 

2  2 

For  a  weak  magnetic  field  B  <4*  nbn>bvo  a  cur¬ 
rent  sheath  may  develop  at  the  beam  surface  to 
exclude  the  external  field.  The  actual  details 
of  the  propagation  for  both  cases  require  a  study 
of  the  evolution  of  the  space  charge  and  current 
sheaths  17] .  However,  an  application  of  the 
virial  theorem  to  the  beam  propagating  in  vacuum 
is  instructive  and  yields  the  global  equation  of 
motion  (6) s 

fd3x  nhm^v.XjL  =  i  fd?x  Xj_.  JxB*  f  d3xP,  (2) 

with  P  =  nbmbv^  +  3pb  +  I(SE)^  +  ( 6B)  2)/4w,  pb 
and  J  =  Jb  +  are  the  beam  pressure  and  net 
current  density  respectively,  j  =  -n  e  vg  is  the 
electron  current  density,  6E  and  «B  are  the  self¬ 
fields.  Note  that  n  -  n.  and  J,  +  J  =0  over 
e  D  be 

the  bulk  of  the  beam  except  in  the  sheaths,  so 
that  the  first  term  on  the  RHS  of  (2)  is  almost 
zero.  Thus  it  is  clear  frcm  (2)  that  the  beam 
must  expand  in  a  vacuun  magnetic  field  regardless 
of  initial  conditions. 

In  what  follows  we  treat  the  injection  and 
propagation  of  ion  beams  into  a  pre- ionized 
plasma. 

Beam  Driven  Return  Current 
The  injection  of  a  beam  of  charged  particles 
into  a  plasma  is  accompanied  ty  a  return  current 
driven  by  the  inductive  electric  field  at  the 
beam  head  (8).  This  current  decays  (9)  on  a  time 
2  2 

scale  T  *  woa  /c  where  o  is  the  conductivity, 
and  a  is  the  beam  radius.  However  if  the  beam  is 
diverging  or  converging,  the  beam  radius,  and 


consequently  the  beam  inductance,  also  varies, 
leading  to  a  source  of  inductive  voltage  that  can 
Irive  or  decelerate  the  return  current.  In  the 
early  stages,  or  equivalently  for  short  distances 
behind  the  bean  head,  we  may  neglect  the  motion 
of  the  background  plasma  ions.  From  Ampere's  and 
Faraday’s  laws  and  employing  Ohm's  law  j  = 

=  -(o/c)  3A/3t,  we  obtain  (10], 


3A 


1  3 


3A 


k.  i  o  on  o  .  ,  .  /  \ 

TrofTFr7r  "3  3b(r'  t_z/vo) ' 


(3) 


where  the  magnetic  field  B  =  VxA  z.  Let  the 
injected  beam  profile  be  given  by: 

=  x  [ I^/wa2 ( t ) 1  exp-(r2/a2)  0  (t-z/vQ)  and 
U(x)  is  the  unit  step  function.  Solving  (3)  by 
standard  techniques  we  obtain: 

jz  =  -|J(t.)  + 
t  2,  2. 

/  dto  <i/a>  ,;TT-/rr  J(r)  0(to_2/von  (4a) 

O 


•  2  2 
with  t*  =  z/v  ,  a  =  da/dt  .  t  =  t-t  .  T  =  no  /c  , 
o  o  o 

I.  exp|-r2/a2(t  )  [1+tAI  l 

J(a(t  ) ,  t-t  )  =  - - - - -  ,  (4b) 

°  °  na^(to)(l  +  tA) 


and  J(t*)  r  J(a(t*),  t-t*). 

The  first  term  jj  of  (4a)  is  the  usual 
response  obtained  when  the  beam  radial  profile  is 
constant,  i.e.,  a  =  0.  The  second  term  jjj 
depends  on  a.  For  an  expanding  beam  with  a  >  0 
this  component  has  the  sign  of  jj  but  for  a  <  0, 
i.e.,  a  converging  beam,  it  apposes  jj.  Thus  in 
a  converging  beam  the  net  return  current  is  dimi¬ 
nished  leading  to  increased  focusing  by  the  self- 
magnetic  field.  The  radial  distribution  of  jjj 
is  different  from  jj  ;  jrl  peaks  near  the  beam 
edge  since  it  is  driven  by  bean  expansion/ 
contraction.  In  the  limit  r  «  T  we  find, 

j  I J  =  -(Ijy/wr2)  x 

2  2  2  2  2  2  2 
f r( 2,  r  /u  tA)  -  r ( 2,  r  v^/z  u  )),  (5) 

where  r(p,c)  =  f  dt  t^*-*  exp(-t)  is  the 

in'  jnpl'-'.u  garrmS  function  and  u  :  a  Is  taken  as 


constant.  The  magnetic  field  is  given  by 

Bs  =  -3A/3r  =  r-1  f  dr  r  (jb  +  jf  +  jn>. 
o 

Beam  Dynamics 

The  momentum  conservation  equation  for  an 
axi-symmetric  beam  (3/36  =  0  =  3/3z)  with  trans¬ 
verse  pressure  P,  density  n^  and  radial  velocity 
u  is  expressed  as: 

It  (nb"\,u)  +  F  7r4rnbnV>u)  =  "  I?  '  F  jbV  (6) 

Integrating  this  equation  over  r  dr  from  0  to  « 
and  recognizing  that  the^number  of  beam  particles 
per  unit  length  =  2*  f  dr  r  is  conserved 
we  obtain, 

Vb  3E  11  -  2"  "  dr  P  -  2c‘2  J  dr  rb  F3F'  (7) 

o  o 

»  r 

with  N.  u  =  /  dr  r  n.u,  I.  =*  2n  f  dr  r  j.  and  the 
b  ^  b  b  b 

net  current  I  =  2*  T  dr  r  (j.+j°);  j  is  the 
6  b  Jp  p 

induced  plasma  return  current.  Let  the  pressure 
2  2 

P  =  nbm5v0( 46 )  i  where  AB,  the  beam  divergence, 
varies  adiabatically  with  exponent  y  =  2,  so  that 

(A©)2  =  (A9>2  a2 /a2.  (8) 

The  radial  scale  length  or  bean  radius  is  a. 
Recognizing  that  the  flux  «  /  dr  B0  is  conserved, 
(7)  may  new  be  written  as,  ° 

NbV'  =  ^bVo7*3  -  (9) 

where  a  =  u ,  I  =  fl.  ,  f  is  the  fractional  current 
h  ——2  *  2 

neutralization  of  the  beam,  TTr/a  =  /  drfl^/r  and 
2  2  b  '  b 

T  =  m.vz(A9)  .  At  equilibrium  £'=  a,  and  with 
O  D  O 

no  neutralization  (f*l)  we  obtain  the  Bennett 
pinch  condition  =  2NbTQ.  0,1  t*le  contrarV  at 
injection,  f  <  1  and  the  Bennett  pinch  condition 
is  not  satisfied.  Choosing  a=0  at  t-0  and  f-0, 
(9)  furnishes  the  solution 

a(t)/a  -  (1  +  t2/t2)1/2  (10) 

o  b 

with  t  »  a  //7  v  ( A# )  .  If  f  *  0  the  first 

DO  0O 

integral  of  (9)  is  given  ty: 


(11) 


\  G2;  *  (a)2  =  < TQ/mb )  (1  -  a2/a2)  - 
(TT^/tJ^ir^c2)  In  a/aQ. 


The  radius  at  wMch  the  self-pinching  force 
finally  catches  up  with  the  pressure  term,  i.e., 
where  a  =  0  is 

af  ”  *0le*P  WX!  =  aQexpn.  (12) 


The  beam  expands  to  radius  a^  in  a  time 
tb  *  (  n/2)  ^//2aoen(tt^/Ubmbc2]  E^2*(n),  where  * 
is  the  error  integral.  If  tb  is  much  less  than  T 
we  are  justified  in  treating  f  «  t/T  as  slcwly 
varyiny  in  (9).  If  the  final  radius  is  not  to  be 
excessively  large  conpared  to  the  injection 
radius  a  ,  it  is  imperative  to  reduce  the  initial 
divergence  (A0>o  and  increase  f  by  re&icing  the 
resistive  decay  time  T  by  inject  ion  into  an 
initially  un-ionized  medium. 

If  the  beam  is  injected  into  a  highly 
resistive  mediim  such  that  initially,  the 
magnetic  force  exceeds  the  pressure  force,  the 
final  self-pinched  radius  af  is  less  than 
and  is  given  by 

(af/aoHn(ao/af)  =  tyr^/FT^.  (12a) 


Finally  one  might  argue  that  since  eventually  all 

tlie  return  current  will  decay  for  a  long  enough 

beam  L  >  v  T,  f  *  1  and  the  tail  end  of  the  beam 
o 

will  eventually  self-pinch.  In  this  situation, 
one  has  to  consider  other  coupe ting  processes 
such  as  kink  instabilities,  etc. 

Plaana  Response  to  Beam  Expansion 
So  far  in  determining  the  return  current  j 
and  its  influenoe  on  beam  dynamics  through  the 
current  neutralization  factor  f  we  have  ignored 
the  response  of  the  background  plasma  ions.  In 
the  same  geometry  of  the  previous  sections  we  now 
include  the  plasma  response  in  the  pressureless 
fluid  approximation.  Thus  with  n^<  n,  the  plasma 

ion  density. 


iz  =  0  (Ez  +  c  W  = 


-£<*?  +  v  i*  ) 

c  St  r  Sr 


o  dA 

=  "  c  at  ' 
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=  0. 


EVr,t)  far  (14) 
+  (i-f)jb)ig,  (is) 

(16) 


Vte  have  three  equations  (12),  (13),  and  (14) 

to  determine  the  vector  potential  A,  the  plasma 

density  n  and  the  plasma  ion  radial  velocity  v^, 

given  j.  (r,t)  =  n,  (r,t)  ev  ,  the  beam  current 
b  b  o 

density  profile. 

Infinite  Plaana  Conductivity  o*» 

This  is  a  particularly  simple  case  for  which 
we  derive  an  exact  nonlinear  solution  of  the  form 
where  the  driving  beam  current  profile  is  jb<r,t) 
=  jjjfr')  with  r1  =  (r2  -  2at)E//2  with  a  con¬ 
stant.  In  other  words  we  regard  the  beam  to 
expand  in  such  a  manner  that  the  local  Lagrangian 
beam  particle  velocity  dr/dt  =  -a/r.  All  quan¬ 
tities  n,  v^  and  A  dept.  only  on  r*  then  since 
ct-k»  and  j  must  be  finite,  iron  (13), 


“  ♦ 
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.  3A 
vr  "Sr" 


(17) 


Note  that  3 A/at  <*  -(a/r' ) 3A/3r’  and  3A/3r  = 

(r/r' )3A/ar' .  The  conclusion  from  (17)  is  that 
either  3A/3r  =  0  or  vf  =  a/r.  If  3A/3r  =  0  then 
B0  =  O 


Jz  +£jb  "  °' 


(18) 


i.e.,  current  neutralization  is  perfect  and  is 

maintained  as  the  beam  expands.  This  is  not  a 

surprising  conclusion  for  o  *  •».  Indeed  if  (18) 

holds  the  beam  may  expand  in  any  arbitrary 

fashion.  On  the  other  hand  if  vf  =  a/r  then  the 

continuity  equation  (16)  is  satisfied  for 

arbitrary  density  profile  n(r').  Furthermore, 

dv  /dt  also  vanishes  and  hence  because  B  *  0 
r  8 


jz(r*>  =  0, 


(19a) 


i.e.,  there  is  no  return  current.  From  (14)  with 
B  ( r * )  *  (r'/r)  B  (r),  we  obtain, 

O  “ 


rB0(r,t)  -  rB9(r,0)  +  ll  /  dr'  r£  f  jb<r’),(19b) 
2  1/2 

where  R  *  (r  +  2at)  .  In  this  solution  the 
plasma  move  radially  outward  with  »  a/r  and 
the  electric  field  vanishes  in  the  frame  moving 
with  this  velocity  and  hence  )  ■  0. 


Finite  Conductivity 

We  new  treat  the  more  urportant  case  of 
finite  conductivity.  In  order  to  sinplify  the 
analysis  of  (13)  to  (16)  we  proceed  to  the  limit¬ 
ing  case  of  slab  geometry  where  r  =  R  +  x 
with  x  <  R.  In  this  limit  (13)  to  (16)  become 

iz  3  -( o/c) ( 3A/3t  +  v^A/Sx),  (20) 


4b 

c 


(fjb  +  V' 


(21) 


mtk  Vr*4v  =  -|?fjb  +  lir4}' 

3x 


3n  .  3 

+  « —  nv  =  0. 
7F  3x  r 


(22) 

(23) 


Let  fjh(x,t)  =  fjb(x‘)  with  x'  =  x-ut,  u  constant 
and  n,  A,  and  v^_  have  the  same  functional  depen¬ 
dence.  From  (20),  (22),  and  (23)  we  obtain. 


j  =  (<j/c)  (u-v  )  3A/3x*,  (24) 

z  r 

n(u  -  v  )  =  n  u,  (25) 

r  o 


v  =  (8wn  mu)~1{(3A/3x*)2+  2  /  dx  j.dA/dx),  (26) 
r  °  x'  b 

with  3A/3x'  =  0  »  v^  at  x'  =  •>.  Subtracting 
these  expressions  in  (21)  we  obtain 


^  =  b  (b2  -  1)  +  fJb  +  2b  /  dx'  Jbb,  (27) 

x 

where  b2  3  B2  /wn  mi2,  B  =  -dA/dx' ,  x  =  x'/S« 

2  0  °  9  2  -1/2 
5  =  c  /4bou  and  J.  =  j.(c/au)(8nn  mu  )  .  In 

D  D  O 

the  region  outside  the  beam  where  J.  =  0,  (27) 

b 

has  the  solution 


.  .l-b.1/2  _  . 
b(I+B>  b 


1-b  , 
_ o.  1/2 

oT+5- 

o 


exp-x 


(28) 


where  b  »  bQ  at  x  *  0  or  x'  *  R  and  b  <  1  every¬ 
where.  Furthermore,  no  solutions  exist  with 

2 

acceptable  boundary  conditions  for  bQ  £  1. 

Indeed  if  bQ  W  then  the  density  n  *  no/(l-b2) 
will  become  infinite  at  some  finite  x»  i.e.,  the 
laminar  solution  ceases  to  exist  because  the 
fluid  particles  are  suffering  a  reflection  at 
this  point  by  virtue  of  the  strong  magnetic 
field. 

Let  f Jb  the  net  beam  current  be  located  in  a 
thin  shell  of  radius  a  and  small  thickness  4  «  a. 
In  this  case  the  solution  within  the  beam  is 
approximately  b  *  0  and  Bq  *  4»f j^A/c.  Then  the 


conclusion  from  the  above  considerations  is  that 
for 

2  1/2 

B  <  (8»n  nu  )  '  ,  (29) 

o  o 

as  the  beam  expands  with  velocity  a  =  u  a  laminar 

disturbance  is  set  up  in  the  surrounding  plasma 

2 

that  typically  extends  for  a  distance  c  /4bou 

beyond  the  expanding  beam.  If  f increases  to  a 

value  such  that  B  =  /8xn  nu2  then  the  distur- 
o  o 

bartce  cannot  be  laminar  and  the  fluid  is 
reflected  ty  the  magnetic  field  of  the  expanding 
beam.  If  this  situation  obtains  then  the  rate  of 
momentum  gain  by  the  plasma  ions  2BanQm(a)2  has 

to  be  balanced  by  the  magnetic  pressure  force 

2  2  2 
(B  /8n)2Ba,  i.e.,  B  =  8nn  mu  ,  which  is 
o  o  o 

precisely  the  condition  (29)  derived  above.  Soon 
after  injection  Bq  is  small  and  u  is  large  so 
that  (29)  is  satisfied  and  a  laminar  disturbance 
propagates  through  the  plasma.  As  Bq  increases 
with  f  and  a  decreases  with  time,  ultimately  (29) 
is  violated.  During  this  last  phase,  the  mag¬ 
netic  pressure  term  in  (9)  should  be  replaced  by 
2*anom(a)2.  If  by  this  stage  the  pressure  term 
is  negligible  then, 

a  3  uQ  exp  -  (a2-a2)/a2,  (30) 

where  a*  3  Nh^/nn^m  and  uq  =  a  at  a  3  aQ. 
Propagation  In  Ionized  Current  Channels 

One  of  the  critical  issues  in  Light  Ion 
Inertial  Fusion  (LIIF)  schemes  is  the  transport 
of  multi -meg  arrpere,  multi-megavolt  light  ion 
beams  over  distances  2-5  m  to  the  target  pellet 
from  the  injection  port.  One  promising  possibil¬ 
ity  involves  the  use  of  preformed  z-pinch  plasma 
colim  channels  to  guide  the  beam  to  target 
[ID  .  After  such  a  channel  is  formed  by  one  of 
several  techniques,  the  azimuthal  magnetic  field 

B.  confines  the  beam  ions  radially.  Space  charge 
8 

and  current  neutralization  are  provided  by  the 
high  plaana  density  in  the  channel 
(n  >  1017cm-3).  The  propagation  of  bean  ions  in 
such  channels  has  been  studied  (i)  using  one 
dimensional  r-dependent  models  [11,121  and  (ii) 
two  dimensional  (r-z)  models  employing  a  hybrid 
2-1/2  developed  simulation  code  CIDER  developed 
at  Cornell  University  f 131 .  These  2-d  simula¬ 
tions  show  that  5-7  MeV,  1-2  MA  proton  beams  with 
initial  divergence  68  of  order  a  few  degrees 
propagate  with  an  energy  efficiency  of  *70  I  in  a 


1 8 

4  m  ionized  hydrogen  channel  of  density  *10 
cm  3  without  revealing  any  gross  axisyimetric 
instabilities.  The  energy  loss  results  from: 

(i)  the  direct  collisions  between  the  beam  ions 
and  the  plasma  ions  and  electrons  and  (ii)  the 
Ohmic  losses  suffered  by  the  return  current  [10]. 

We  will  consider  here  the  dynamics  of  beam 
propagation  in  such  channels.  Prior  to  the 
injection  of  the  beam  into  the  channel  it  is 
assumed  that  an  equilibrium  exists  between  the 
magnetic  and  pressure  forces,  i.e., 


B 

3P  =  -  9  »  r  o 

Ir  47F3rrB0' 


(31) 


For  high  enough  channel  conductivity  the  flux 
t|>  =  f  drB  is  conserved  during  the  expansion  of 
channel .  Thus 

a"=  K/a  :  K  =  ( 1^/31^) (B2/4nrm) .  (38) 

The  solution  of  (38)  is  given  by, 

a  =  f2K[lna/ao)  11/2,  (39) 

(a/aQ)  F[ (In  a/ao)1/2]  =  (K/2a2)  t,  (40) 
2  x  2 

where  F(x)  =  exp(-x  )  /  dy  exp(y  ).  For 

a  =  a  +  4a,  5a  <  a  we  obtain 
o  o 


2 

and  I  =  na  j  is  the  channel  z-current. 
c  c 

Furthermore,  the  beam  must  be  injected  before  any 
gross  instability  of  the  channel,  viz.  the  kink, 

sets  in  on  a  time  scale  t.  =  y  *  where 
1/2  k 

Y  *  (Bo/(4nnm)1/z  ]k  IndAa);  Bq  =  2»jca/c,  a 
is  the  channel  radius,  k  =  2n/X  and  y  is  the  kink 
wavelength.  The  injected  beam  is  expected  to 
satisfy  the  equilibrium  relation 

5  -  -  £  jbBe  +  W  (32) 

where  E  is  the  Hall  radial  field  required  to  set 

r  2  2 
up  the  return  current,  and  P  *  .  If 

the  injected  beam  does  not  satisfy  (32)  then  it 

may  either  expand  or  contract  until  it  settles  to 

an  equilibrium  radius  determined  by  an  equation 

similar  to  (12)  with  replaced  by  I^^a/2. 

The  channel  plasma  dynamics  is  governed  by: 

im  ^  =  -Vp  +  I  jc  x  B  +  ^  jpX  B  -nfee  E,  (33) 

j  =  <j(  E  +  ivxB  -  j  X  B  /nec ) ,  (34) 

p  c  p  0 

cVx  B  ■  4x<  ^  +  jp),  (35) 

ne  =  (36) 

Inmediately  after  beam  Injection  the  induced 

plasma  current  j^  *  -  j^*  -n^ev^.  From  the  r- 

conponent  of  (34)  and  j  =  0  we  obtain  the  Hall 

electric  field  E  *  (n. /n„)  v  B./c.  Thus  for 
r  d  0  o  “ 

nb<  n  ,  E^  may  be  neglected  in  (33).  Integrating 
(33)  over  rdr  and  using  the  continuity  equation 
we  obtain 

m  3F  \  1  ""  *  ’  dr  W*  (37) 


6a/ a  =  t  /t  ?  (41a) 

o  c 

t  =  (31  /lM2)a  /v  ,  v2  =  B2/4wn  m.  (41b) 

C  CDOAAOO 

Now  since  the  beam  ion  betatron  time 

t  =  a  /v  A9  C  t  ,  the  beam  contained  by  the 
B  o  o  c 

channel  magnetic  field  expands  adiabatically  with 

the  channel  radius  but  satisfying  (32).  From 

(40)  for  a/a  =  2,  t/t  =  0.9  and  for  a/a  >  1 
o  c  o 

the  beam  and  channel  expand  almost  linearly  with 

time.  If  the  beam  expansion  is  to  be  limited 

during  propagation  the  beam  pulse  width  r  <  t^. 

The  channel  expansion  may  ultimately  be  solved  by 

the  dense  neutral  gas  surrounding  the  channel 

that  absorbs  the  momentum  through  ion  neutral 

collisions  provided  v.  t  >  1;  v.  is  the  ion 
in  c  -  in 

neutral  collision  frequency. 

Kink  Instability 

After  a  time  t  the  return  current  decays  a 
little  and  the  net  beam  current 
I  =  flb  (f*tc/T,  T  *  »aa2/c2)  Increases  the 
magnitude  of  the  magnetic  field  causing  the  beam 
to  pinch.  This  could  lead  to  non-axisymnetric 
instabilities  not  covered  by  previous  1-d  and  2-d 
simulation  studies  [11,12,13).  Since  the 
resistive  hose  grows  on  the  time  scale  T  and  the 
pulse  width  r  <  T,  we  need  not  be  concerned  with 
this  particular  instability.  On  the  other  hand, 
the  kink  instability  occurs  also  in  the  infinite 
conductivity  limit  and  we  must,  therefore, 
examine  this  case  next. 

For  the  purpose  of  this  calculation,  we 
shall  assume  that  either  the  channel  has  stopped 
expanding  or  that  the  expansion  rate  is  much 
slower  than  the  growth  rate  y  if  any,  of  the  kink 
instability,  i.e.,  y  >  a/a  •  t~3.  bet  the  plasma 
displacement  be  described  by  C*exp  i(kz+9-wt) . 


Then  from  the  perturbed  mhd  equation  for  plasna 
motion  we  obtain 

2  3  -  ,2  *  *  3 

— w  fd  xnm  Id  =  -  e/fn^E  •  SB  +  E  *B  )d  x  - 

f  d3x  E  •  78p  +  i  /  d3x  E  *jp  x  «B  + 

(42) 

^  /d3x  E  •  (V  x  «B  )x  B  -  i  / d3x  E  •«  jb  x  B  . 

How  with  S£  respresenting  the  perturbed  beam 
distribution  we  write, 

1  3  *  3  * 

Ij  =  -  —  fd  x  E  *5  jbx  B  -e/d  xSn^  E  •  E 

33*  1 

=  -e  f  d  xd  v  E  •(£•»•  -yx  B  )4f 

c 

=  -  mp  f  d3xd3v(  E  •  ^)8f 

=mb/d3xd3vf <  E  "  v  >3t  5f+  <  v  ‘  3t  E  >4f )» 
(by  virtue  of  Liouvilles  theorem) 

=  fd3xd3vfe(  E  *v  )(  «E  +  i  v  x  «B  )• 

+  (  v  *  3E  5  )sf) • 

(enp laying  Vlasov  equation) 

=efd3xd3v(nbE  •  «E  +  E  •  jfax  «B  )  + 

/  d3x«f(  v  •  ^  E  )• 

Employing  the  above  expression  in  142]  we  obtain 
[14] 

m2  fd3xrm|  E  |2  *  -/d3x8p(7«  E  ) 

-  i  f  d3x  E*-  (  jp+  jb)x  SB 

+  ^  /d3x|  SB  |2  +mbfd3xd3vSf (  v  •  ^  (*).  (43) 

* 

For  the  kink  stability  we  set  7  •  E  =0  thereby 
annihilating  the  first  term  on  the  RHS  of  (43). 
The  last  term  of  the  RHS  can  be  conputed  by 
solving  for  8f  from  the  Vlasov  equations 

of  «  —  (8A  -  -  l(  to  —  + 

3 vz  3e 

3f  fc 

k  |£_)  /  df  vz*Az|.  (44) 


(Where  8A„  =  EB„ ,  E  =  E  (r  +  i9)  expi 
(kz  +  9  -  ut),  e  =  v  /2,  e  is  constant.) 
Substituting  8f  in  I ^  and  enploying  betatron 
motion  of  the  beam  ions  to  evaluate  the  time 
history  integral  we  obtain, 

,32  * 

Ij  =  nib  (  d  xd  vSf  v  •  (-ito  +  v  .  7)  e 


IKB,(kv  -to)2 
b  o 

- 5 - 7  • 

(kv  -toP  - 

O  0 


In  arriving  at  this  expression  we  have  chosen  the 
beam  distribution  f  -  exp-fmhAhHe-v^v^)  and 
B  =  B*r/a  is  determined  by  the  net  current 

v  * 

density  =  jp  +  jb  =  z(jc  +  fjb>,  where  f-t^. 

The  betatron  frequency  of  the  beam  ions  is 
1/2 

Ug  »  (eB^v^^ADhac)  .  The  remaining  terms  of 
(43)  represent  the  nhd  contribution  from  a  psj 
channel  carrying  the  net  current  Ic  +  flp  and  can 
be  treated  in  the  conventional  manner.  Finally 
we  obtain: 


,  B*  ,  ,  ,  I.  (kv -a2) 

“  * - 3f“k  a  inj^  +  < ;  — ) - 2 - 3-t  • 

4wma  *c  rlb  (kvQ-(i>  )-&i  8 

(45) 

where  B*=2(I  +fl.  )/ac.  The  beam  term  is  stabil- 

CD 

izing  for  perturbations  with  wavelengths  satis¬ 
fying  the  inequality  [16] 


.2  ^  2.  2  _  .2 

k  >  Vv0  5  V 


(46) 


For  stable  propagation  we  also  require 


2  2 


kaln(lAa)  <  ipA^  +  flb>- 


(47) 


Identifying  the  longest  unstable  wavelength  with 
the  beam  pulse  length  L  -  vqt,  stable  propagation 
is  achieved  if 


(2wa)2(Ic+  flb)/lb  <  v2x2  <  (2wvo/toB)2.  (48) 

This  result  is  acocrpanied  by  a  caveat.  For 

unstable  long  pulse  lengths  the  growth  rate  y 

from  (45)  is  approximately  yt  -  2w/«8t*s  t* 

is  obtained  from  (41)  with  I  +  fl  replacing  I  . 

c  b  c 

The  basic  justification  for  the  analysis  that 
requires  yt  *  >  1  may  not  hold.  However,  the 
treatment  of  the  stability  of  an  expanding 
current  column  is  nontrivial  but  it  is  unlikely 
that  an  exact  treatment  will  lead  to  a  more 
pessimistic  picture  because  an  esqoanding  channel 
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reduces  the  free  energy  available  for  kinking. 

In  any  case,  the  nunber  of  e-folds  yt  -  t  /t  * 
computed  from  the  above  analysis  is  less  than 
unity  in  a  dense  channel, 

Filamentation  Instability 

Although  generally  speaking  the  beam  pulse 
time  i  «T  the  return  current  resistive  decay 
time,  this  may  not  be  case  for  a  filament  of  the 

bean  of  width  A  «  a  for  which  the  decay  time 
2  2 

Tfi  -  *oA  /c  .  The  decay  of  the  return  current  in 

the  filament  is  acconpanied  by  a  magnetic  field 
1/2 

B  =  2*enbv  it/T)  A/c  at  the  filament  surface. 
The  filament  will  self-pinch  via  the  filamenta¬ 
tion  instability  (17, ]R),  at  say  t  =  t/2  if 

n  =  8»P/B2  =  4nb(A9)2o/nbe2T  <1.  (49) 

Note  that  n  is  independent  of  A  and  is  identical 
to  NT  c2/FT?  of  (12)  for  f  =  t/2T.  If  the  beam 

O  D  T  — 

satisfied  (32)  at  injection  then  ( A 0 )  ~ eBoa/n^cvo 
and  n  =  4  I /f^I^.  If  n  >  1  the  beam  will  prop- 
pagate  without  filamentation.  On  the  other  hand 
if  n  <  1,  filamentation  will  proceed  until  the 
bean  divergence  increases  through  nonlinear 
evolution  such  that  the  state  n  =  1  is  reached. 
Thus  the  filamentation  instability  does  not 
appear  crucial  to  propagation  in  Ionized  current 
channels. 

In  conclusion  it  appears  safe  to  predict 
that  intense  ion  beams  of  suitably  short  duration 
may  safely  be  propagated  in  straight,  well  formed 
current  carrying  channels  of  sufficient  mass 
density.  It  is  amusing  to  note  that  ultimately 
the  confinement  mechanism  for  the  beam  relies  on 
channel  plasma  inertia. 
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